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ON A CLASS OF BILATERAL GENERATING FUNCTIONS 
FOR HERMITE POLYNOMIALS OF TINO VARIABLES 


G. K. PANJA* AND D. K. BASU** 

(Received 20.11.2000 and in revised form 10.11.2001) 

ABSTRACT 


A new class of bilateral generating functions for Hermite polynomials of two variables (x, y) 
is obtained. Applications of our result are pointed out. 

Key word & phrases; Special functions & bilateral generating function. 

Mathematics Subject Classiflcation (2000): .'13C45 


INTRODUCTION 

Noticing the existence of the following type of generating function of Hermite 
polynomials of two variables (a:, y) [5] 


1-^1 w „ 




n\ 


( 1 ) 


where a is a nonnegative integer, we are led to investigate a more general class of 
generating function by Lie group - theoretic method. 

The main result of our investigation is the following theorem ; 

Theorem : If there exists a liner generating function of the form 


n^O 

then the following new bilateral generating function [3] will exist. 
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ON A CLASS OF BILATERAL GENERATING FUNCTIONS 


GROUP THEORETIC DISCUSSION 

For // (x, y) we have the following differential operator [5] 


' d d 2yz d 

a = 2 y — - - -+2jt 

X dy X dz 


such that [3] 


B H„,,{x,y)z'’*' 

The extended form of the transformation group generated by B is [3] 


(e.phB)f(x.y.z) = ^xp 241 -^ / 

1 2x / V X 




DERIVATION OF THE GENERATING FUNCTION 

Let us consider the generating function 

G{x,y,t) =f^a„H,Jx.y)t’‘ 


Multiplying both sides by , we get 


r Gix,y. t) 

/i=0 

Replacing / by tz, we get 
rz‘^G(x,y,rz) = £ 

nsBf) 

Operating both members of eqn. (9) by exp (bB) 

{exp(f»R)}[r“z“GU,y.rz)] = {exp(f;B))£fl,//„,,(x.y)z“"V^ 

/i*0 

The left hand side of the last equation becomes 
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byz. 


exp 2^X20-:^) 
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On the other hand using (5) the right hand side is reduced to 

t=0n=0 

Hence we get 




k=0 n=0 


( 10 ) 


Now putting ^) = 1 we get eqn. (3) using Lemma 10 of Rainville [6, pp. 56, 57] 

APPLICATIONS 

1 ^a+n(^*y)^ 

If we consider = —: then G (x, v, 0 = S-- 

rt! ,.o n' 

Also from [5] we have 
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Again, we obtain from [5] 
G 


41-—)’',>'(l--)‘".tz(l--)' 


yz(\ 


_ ^2xzi-2yK^ -yi^z^ 
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From equations (3) and (12) we obtain 

exp [2x2(1+ + + H- 


( 12 ) 
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Now applying the relation HJ^x,y) = H„ (x / yfy) we get on 

putting y = 1 

exp [2xzil + t)-z\l + tf]H^ [X - z(\ + 01 = S Z (14) 

which is a new bilateral generating function for Hermite polynomial of a single 

variable. 
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BHASKARA'S methods for computing 

MEAN PLANETARY POSITIONS^ 

S. L. SINGH* AND RAMESH CHAND** 

(Received 15.01.2001) 


ABSTRACT 

In this paper we discuss several methods from Siddharua ^iromarii (SiS) 
for computing mean planetary position of planents. A few applications are 
also discussed. 


INTRODUCTION 

In his Grahaganita (c. 1150 A.D.), Acdrya BKdskara (generally called 
Bhdskara II) ( b. 1114 A.D.) (Tienceforth called Acdryain this paper) has 
given sevreal computational methods for the determination of planetary 
position at a particular place and time. The classical presentation of 
Acdrya's highly respected formulae are apparently not very friendly to a 
modern astronomer. In this paper, we attempt a systematic presentation 
of these methods which can be used very conveniently even by a 
sophomore in astronomy. Indeed, here we present six methods concerning 
the mean position and illustrate them by examples. 


t A part of the talk given at the National Symposium on Ancient Science in India, organized 
during the third annual conference of the Vijnuna Parishad of India [dedicated to the 800th death 
anniversary of Bhaskaracarya (1114-1193 A.D.)], H.N.B Garhwal University, Srinagar (Garhwal), May 
25-26, 1993. 

* Department of Mathematics, Gurukula Kangri University, Hardwar-249 404. 

♦♦ Department of Mathematics, V. M. I. C. Sec -5, B.H.E.L., Hardwar-249 403. 
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BHASKARA'S METHODS FOR COMPUTING MEAN PLANETARY POSITIONS 


COMPUTATION OF MEAN POSITION 


1. First Method 

VRIwnl^ 11F:I 
^!nf$iy[;>iR ^ MUW: 11 


This means (see [ 2, p. 31]) 


[SiS, ] Grahanayanadhyaya ] 


The ahargana multiplied by the number of sidereal revolutions of the planet 
and divided by the number of civil days in a kalpa gives the mean position of the 
planet at the time when the mean sun is very nearly at the eastern horizon in Shrilanka. 

This, using notations, may be presented as follows ; 

Mean longitude of a planet in Shrilanka = (R X A ) -^ C, wherein R is the 
sidereal revolutions of a planet in a Kalpa. A = ahargana and C = civil days in the 
present Kalpa. For details of ahargaija, Kalpa and other technical terms one may 
refer to [1-11] and [13], 

Example. Compute the mean longitude of the planets on May 25, 1993. 


Following standard techniques of calculating ahargana (cf. Lahiri [7] and also [13]). 
A {ahargana) on May 25,1993, A=5168 X 360 + 3 X 30 + 30 + 25 = 1860667 civil 
days. Therefore, on May 25,1993, mean longilude of the Sun. 

= (432 X 10' X 1860667) / 1577916450000 revolutions 

= 5094. 110933 revolutions 

= 5094 revolutions /1 Rds'i / 9°/ 56V 9.68 
mean longitude of the Moon 

= (577533 X 10* X 1860667 )-i-1577916450000 revolutions 
68102.24896 revolutions 

= 68102 revolutions / IRdsi / 297 37'/ 32.16 
mean longitude of Mars 


229828522x1860667 

1577916450000 


Dhruvakas 


= 271.0120338 Revolutions +11 Rasi / 29®/3'/50" 
= 272 Revolutions / o Rdsi/3°/ 23V 45.4048". 
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Similar computations may be made for other planets (see, for example, [13]). 
2. Second Method 

131490000000 

csranmi 


[SiS, Grahcawyanadhyaya] 
The avama-k^a divided by 13149 X 10’ in (degrees) is to be added to twelve 
times the elapsed tithis and the result added to the sun's position gives the position 
of the moon. Conversely the position of the sun can be obtained from the position of 
the moon (see,[2, p. 32]) 

This means Moon's longitude = (x+12 Y)'’+R“ /13149 x 10’ 

(wherin x = sun's longitude, Y = elapsed tithis up to the concerned day, 

R = avama s'easa). 

For the computation of elapsed tithis till May 25, 1993, we as fallows: 

D (elapsed solar days till May 25,1993) 

= {[3179 + (1993 - 78)] x 12 + 2 } x 30 + 3. 

= 1833903. 

Elapsed intercalary months 

= DA/S, (where A = adhikatndsaof the Kalpa, S= number of solar days in Kalpa) 
- (1833903 X 15933 x 10') h- 15552 x 10» 

= 1878.83791 

~ 1878 (Discarding the remainder as per rule). 

Y (elapsed tithis) 

ss D+elapsed intercalary months x 30 
s 1833903+ 1878x30 
= 1890243 

Since /? = (Y x ot) T (where Ot = Omitted tithis in kalpa and T tithis in Kalpa) 
we have: 

R = (1890243 X 2508255 x 10^) ^ 1602999 x 10* 

= 29577. 13294 
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The mean longituded of the Moon on May 25,1993 
= (x+ 12Y)« + R‘>/13149 X 10’ 

= (5094. 110933 + 12 x 1890243 / 360) revolutions + (29577. 13294/ 1349 x 

10 ’)« 

= 68102 Revolutions/ 2 RajJi 1575675.2808067". 

3. Third Method: 

^ T?af 

crar erar «fci^ ^svrarw; 11 

[Sis Grahanayahadhyaya] 
Even as the sums or differences of two or more of the numbers of adhimasas, 
avama ( ksaydhds ), lunations etc. give the number of sidereal revolutions of the 
planets, the sums or differences of two or more of the positions of the imaginary 
planets which go by the names adhimasa- graha , ksayahagraha etc. computed out 

of the numbers of these adhirmsas ksayahas etc. give the respective planetary 
positions. (See [2, p. 39]). 

According to this, longitude of the moon on May 25, 1993 
= 13 X Revolutions of the sun + Intercalary- months 
= 13 X 5094. 11093 + 1878. 830791 
= 68102.272911 Revolutions 
= 68102 Revolutions / 3 Ras'i / 8714752.656". 

4. Fourth method: 

9921 

3 5419 cRWni 
11859 

’jErai 'ei 

Iclrf^ gr TF^ %yr|u rif ^d>l : W:l 

[Sis, Madhyadhikara Madhyadhikara ] 
The ahargana multiplied by 11859 decreased by the quotient obtained by 
dividing the product of the ahargana and 9921 by 35419 given the distance covered 
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by a planet in yojanas. These yojanas divided by the circumference of the planet's 
orbit gives the fraction of a revolution and the integral number of revolutions made. 
(See [2. p. 50-51]). 

This means distance covered by a planet in Yojnas 

= Ax 11859-Ax9921 / 35419, 

where A is the required ahargana on the concerned day on May 25,1993 and 
A = 1860667 days. 

Thus the distance covered by a planet 

= 1860667 (11859- 9921/ 35419) Yojanas 
= 1860667 (11858.7199) Yojanas 
The mean longitude of the sun on May 25, 1993 

= (Distance covered by the planet in the required ahargana) 

■r (circumference of the planet's orbit) 

= (1860667 X 11858.7199)- 43314975 Revolutions 

(wherein 4331497.5 Yojanas is the circumference of the Sun's orbit) 

= 5094.110924 Revolutions 

=5094.110924 Revolutions/1 Ras'i/ 9755'/57.504" 

The mean longitude of the moon on May 25, 1993 

= (1860667 X 11858.7199)/324000 Revolutions 

(Wherin 324000 Yojanas is the circumference of the moon's orbit) 

= 68102.24931 Revolutions 
= 68102 Revolutions / 2 Ras'i / 29745'/5.76" 

5. Fifth method 

(V.l) Longitude of the Sun 

f^vjrsFf^ccTerrf^rcrt 

[Sis, Pratyabada s'uddhi mmadhyddhikdra] 
The number of the days in the ahargana is to be diminished by 1 / 60th part of 
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itself to obtain the number of degrees and fractions thereof; then the ahargana 
multiplied by three and divided by 22 gives the minutes and fraction there of. Adding 
the two results, we get the reqired position of the Sun. (See[2.p.64]) 

Accordingly the mean longitude of the Sun 

= {Ahargana - Aharganat6QP) + (3 x Ahargana/22) 

= Ahargana [ (59/60)“ + (3/22)' ] 

The mean longitude of the sun of May 25, 1993 
= 1860667 [(59/60)“ + (3/22)’] 

=5094.124033 Revolution 

=5094 Revolution/ 1 Ras'i//14“/39'/6.768" 

(V.2) Longitude of the Mars : 






[Sis, Pratyabada s'uddhi in 
Madhyadhika ra ] 

" Half of the number of days gives the degrees; half the number of days 
multiplied by 3 gives the number of minutes "; from this subtracting 17th part of 
Ahargana gives the minutes and adding the Drhuvkd at the beginning of the solar 
year, gives the mean position of Mars (See [p.65]). 

Thus the longitude of the Mars 

= (Aharagaria / 2)“ + (3 Aharagar^a / 2 )'- (Aharagar^a 117)' + Dharuvaka. 
The longitude of Mars of May 25,1993 

= (1860667/2 + 3 X 1860667/2 x 60-1860667/17 x 60) + Dhruvaka 
= 2708 Revolutions / 145“/59'/23.8992’'+l 1 Ras'i /29“/3'/50’' 

=2709 Revolutions/4 Ras'i / 24“/62'/73.8992". 
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(V.3) Longitude of Mercury's S'ighrocca : 

^ ^ dd l ljdl 
^trsIH 11 

[Sis, Pratyabada s'uddhi in Madhyadhikdra] 
The Ahargana multiplied by 4 and 4/130, gives respectively the positive degrees 
and minutes to be added to the Dhruvakds of Mercury's s'ighrocca. 

Thus the longitude of Mercury's s'ighrocca 

=(4 Ahargartdf + (4 Ahargana 3/130)' + Dhruvaka 

The mean longitude of Mercury's s'ighrocca on May 25, 1993 

= (4 X 1860667 + 4 x 1860667 x 3/ 130)" + 11 Rds'il2T-l2A'l29" 

= 21151. 172 Revolutions + 11 Rds'H 27“/24'/29" 

= 21152 Revolutions / 58"/2'/0’' 

(V.4) Logitude of Jupiter: 

[Sis, Prafyabada s'uddhi in Madhyadhikdra] 
The ahargana divided by 12 and by 71 gives respectively the positive degrees 
and negative minutes to be added to the Dhruvaka of Jupiter. (See[2, p. 66]). 

Thus the longitude of Jupiter 
= (Ahargana/\2)° - (Ahargar^a/l V)' + Dhruvakds 
The mean longitude of Jupiter on May 25, 1993 

= (1860667/12)" - (1860667/71)' + 11 Kds'H 29"/27'/36" 

= 430 Revolutions / 5 Paj'i728"/16'/1.30244" 

(V.5) Longitude of Venus S'ighrocca: 



^9^ %cf:ll 

[SiS. Pratyabada S'uddhi in Madhyddhikdra] 
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The Ahargana multiplied by 10 and divided by 6 and 155 respectively gives 
degrees positive and negative to be added to the Dhruvaka of Venus to give his 
position. (See, [2, p. 66]) 

Thus the mean longitude of Venus S'ighrocca 
= (10 Aharganaf 6) (10 Aharganal 155)° + Dhruvaka 

- The mean longitude of Venus S'ighrocca on May 25, 1993 
= 1860667 [10/6- 10/155]° + 11 /?aj'i/29°/42V14" 

= 8280. 7462 Revolutions + 11 /?a5'//29°/42714’’ 

= 8281 Revolutions / 8 Ra5'i/27°/20’/17.7536" 

(V.6) Longitude of Saturn 

[SiS,Prafya/?ac(a S'uddhi in Madhyddhikara] 
This means that the longitude of Saturn 

= (2 X Ahargana)' + (2 x Ahargana/S)" + Dhruvaka 
The mean longitude of Saturn on May 25, 1993 

= (2 X 1860667)' +(2x1860667 / 5)" +11 Ras'i /28°/46’/34" 

= 172 Revolutions/308°/58V26.78736"+ll Ras'i +28°/46y34" 

= 173 Revolutions/10 Ras'i/7°/45V0.78736" 

(V.7) Longitude of Moon's Apogee : 

left 

[SiS, Pra/yabada S'uddhi in Madhyddhikara] 

The Ahargana divided successively by 10 and 88 and added gives the degrees 
to be added to the Dhruvaka of the apogee of the moon to obtain its position. 
(See [2, p.67] 

This means that the longitude of Moon's apogee 
= (Ahargana / 10)° + (Ahargarfal 88)° + Dhruvaka 
This means longitude of Moon's apogee on May 25, 1993 
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= 1860667 (1/ 10 + 1/ 88)® + Dhruvaka 
=575.58512 Revolutions + 4 Ras'H5QPl29'l46" 

= 575 Revolutions / 11 /?as'i76®/8721.52" 

(V.8) Longitude of Lunar Node: 

W»RIT:I 

^ ^ TT^II 

[Sis. Pratyabada S'uddhi in Madhyddhikara] 
TheAhargam multiplied by 30 and divided by 566 gives the number of degrees 
to be added to Dhruvaka to give the position of Rahu (lunar node). (See, [2, R68]). 
This means that the longitude of lunar node 

= (30 Ahargaria/566f + Dhruvaka 
The longitude of the lunar node on May 25, 1993 
= (30 X 1860667/566)® +5 Rds'il 3°/ 12y58" 

= 273.94979 Revolutions + 5 Ras'i / 3® 127 58" 

= 273 Revolutions/ 341® / 557 27.84" + 5 Ras'i / 3®/ 127 58" 

= 274 Revolutions / 4 Ras'i /15°/ 87 25.84" 
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AN IDEALISED MODEL OF AN M.H.D. SPHERICAL 
SHOCK WAVE IN SELF-GRAVITATING GAS 

A. K. DWIVEDI 

(Received 20.01.2001 and in revised form 15.11.2001) 

ABSTRACT 

An analytical solution of a shock, in a spherically symmetric model incorporation a self-consistent 
azimuthal magnetic field has been studied using chester-chisnell-whitham (C.C.W.) method. The 
shock is assumed to advance into a self-gravitating and conducting gas of spatially decreasing density 
and pervaded by an idealised spatially decreasing magnetic field. 

Numerical solutions are obtained for shock velocity and shock strength. 

KEY WORDS: Spherically symmetric model, Chester-Chisnell- Whitham method, shock 
waves. 

INTRODUCTION 

Summers [10] has constructed an idealised model of a magnetohydrodynamic 
spherical blast wave applied to a flare produced shock in the solar wind. He has 
omitted gravity in his problem and solved it using the similarity method. Recently 
Singh and Mishra [7,5], Singh and Pandey [8,9], Yadav [12], Kumar and Prakash 
[3] and Kumar and Saxena [4], have discussed different types of problems of 
cylindrical shock waves allowing the effects of gravitational force, magnetic field 
and also in rotating atmosphere. They all used Chester [1]- Chisnell [2]-Whitham 
[11] (C.C.W.) method for their solution. None has attempted these problems of shock 
wave in spherical symmetry in magnetogasdynamics using the Whitham's rule. We 
have investigated the analytical solutions for the propagation of spherical shock 
wave in selfgravitating gas of spatially decreasing density in the presence of spatially 
decreasing magnetic field. The problem has been studied for weak shock with weak 
and strong magnetic field as well as strong shock with strong magnetic field. 


* Department of Physics, R.S.K.D. Degree College, Jaunpur-222(X)1 
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Analytical expressions have been obtained for shock velocity and shock strength 
with other variables also. 

According to Summers [10] and Rosenau and Frankenthal [6] we have taken 

the density distribution Po ahead of the shock to vary as an inverse power or radial 
distance from the point of explosion, i.e. 

p.=p;-" (1) 

where p„ and w are constants. 

The magnetic field //<, is replaced by a physically unrealisable idealised field 
such that lines of force lie on hemisphere whose centre is the explosion point as 
shown in the figure one given in the paper of Summers [10]. The magnetic field 
distribution is assumed to be of similar form, i.e. 




5 


( 2 ) 


and directed tangential to the advancing shock front. Equation (2) is introduced 
purely to give the mathematics of the problem tractable and its in correct form will 
not substantially affect the overall validity of the present model. 


BASIC EQATIONS AND BOUNDARY CONDITIONS 

The equations governing the spherically symmetrical flow of a gas under the 
influence of its own gravitation in the presence of azimuthal magnetic field following 
Summer's, are 


du du I dp ^GM dH ^ p/f ^ ^ ^ 


’ -f- w- 


dt dr p dr r 


dp dp 
—+ M—-i-p 
dt dr ' 


^ r ; 


dp dp 2 
—+u—-a^ 
dt dr 




dt dr 


p dr 
= 0 

= 0 


P'- 


(3) 

(4) 

(5) 


dm 


-47ir^p = 0 


( 6 ) 
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dH dH „ 
dt dr 


(du u 


+ 

dr r 


= 0 


(7) 


where m {r, t), u (r, /), p (r, t) and H (r, t) denote, respectively, the mass inside 
a sphere of radius r, the velocity, the pressure, the density and transverse magnetic 
field at a distance r at time /. p is magnetic permeability and 




The magnetohydrodynamic conditions can be written in terms of a single 

parameter N = — as 
Po 


p, = Np,„H, = NH„u, = {l-l/N)U 
2N 


= 


Pi = Po + 


(Y + l)-(y-l)Af 
2Po(A^-l) 


[a^o+W-y)N + y}] 


( 8 ) 


(y+ 1 )-(y-1)A^ 


2 . 7 ~ * U/x, 1n2 


Uo + 


■biN-\y 


Where 0, and I, respectively stand for the states just ahead and just behind the 
shock front. U is shock velocity, is the sound speed (yPo ! Po)'^and is the 
Alfven speed {\x.Hl / p)'^^. 

1. Weak Shocks 

For very weak shock we take the parameter as 

^ = N = l + e (9) 

KO 

Where e «1. Now consider the two cases of weak and strong magnetic fields. 
Case -1 

For weak magnetic field b] I al«\ under this condition and (9), the boundary 
conditions (8) for very weak shocks reduces to 
p, = Pod + e), H, = Wod + e),u, = eao 


(10) 
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= (l + ^^^8)ao,p, = Po{\ + ye) 

4 

Case - II 

For strong magnetic field bit c^» 1, using this conditions and equation (9), 
the boundary condition (8) reduces to 

u = ^\+^^o. Pi=Po(i+y^) ( 11 ) 

2. Strong Shocks 

In the limiting case of a strong shocks p, / is large and in the presence of 
the magnetic case this can be brought about in two ways. 


Case-I 


Y + 1 

The purely non-magnetic way when A^-> -— 

y-l 

Case-II 

When bo »al or when \i.Hl » Po.that is when the ambient magnetic 
pressure is large compared with the ambient fluid pressure. In terms of Wthe boundary 
conditions (8) now become 

p, = yvp„.yv, = NH,u = {\-\i N)u, 

PL=:X{N)^+YiN) ( 12 ) 

A) fli. 

Where 


XiN) = 


y{y-l){N-lf- 
2N[2-y]N + y] 


Y(N)^ 


N{y + l)-{y-\) 
{y + \)-(y-\)N 


For diverging shocks the characteristic form of system of equations (3) to (7)- 
that is, the form in which equation contains derivatives in only one direction in (r, t) 
plane is 
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j t/jo j Ipc^udr \jLH^(u-c)dr pcGmdr _ 
dp + \iHdH + pcdu— - - —^ — +- - - = 0 


u + c r 


u+c r 


Where 

P P 


ANALYTICAL SOLUTION 

For the final step, substitute the stock conditions (10) or (11) or (12) into 

relation (13). A first-order differential equation in e (r) or is obtained which 
determines the shock. For weak shocks in the presence of a weak transverse magnetic 
field on substituting the stock conditions (10) in (13), neglecting the second and 
higher order terms of e (since 8 is very much less than unity), we find that. 

(dp. da,, 2dr 

2 + - - d 8 + — !2- + -+-i-E 

I ypi) ) . Pi) r yPi) 


ndp ^ d{\xHl) ^ Hi dr ^ Gmdr 
«oVPo 2po Po r r- , 
now in view of hydrostatic equilibrium prevalent in front of the shock 

\r ) Po P,, P„ r 

Equating (14) due to (15) reduces to 

i/:fi t‘ffoY‘<Po 


2YPoA Po 


dCn 2dr 
+ — 2 . +- 


Also we have 

m =- 

3-iv 

27tGp* 2(1 -h) 1 

Inserting 5 = w -1 in the value of p. we have 
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/»o = 


2nGpl 1 (2 - w) 2 
[(w-3)(l-w)'^2(w-l)^' 




Therefore 

and 




(18) 

(19) 


Where 


^,= 


2nGel ,1 (^- 
(h.- 3)(I-«-) 2(»-ir ' 


and 


^2 ~ ^ Pc 


Positivity and finiteness of the equilibrium pressure as defined by equation 
(18) requires that the constant w should obey the inequality 

1 < w < 3 (20) 

Substituting dp„ / p„ and da,,/a„ in equation (16) we get 


— =-[P^(2 + 3w-) + lOw - 20]— 

E 8 r 

Integrating the above equation we get 
e = kr[p-(2 + 7tw)+\0w - 20]/8 

2 

where P = ~i~ and k is constant of integration. 

y^i 

With the aid of equation (10), we can write 

U ^(2-h)/2 y + 1 ^^[p^2*3H)t6H^-l2]/8 


( 21 ) 

( 22 ) 


(23) 
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U , Y + 1 , fp‘(2 + 3K)tlO»-2ul/ 

-- 1 + -i- kr'- ■' 

a„ 4 


Similarly, for a strong magnetic field by use of shock condition (11) in (13) 


we get 

- 1/2(1 I I I 

e 2MWr r ^//o 


Substituting the value of dp^, dHg and db„/b^^ in this equation and integrating 


we get 


[(|/P’)]( 2 + 3 m>)+ 10 u- 20)/8 


e = kr 

Where k is constant of integration. 

With the aid of equation (11) we can write 


^ - ft 


k. 


|(l/(j‘)(2tlH)+«H-l2]/8 


— = b l + (3/4)j 

«() L 


[(1/P*)(2 + 3h)+IOh-2()]/K 


For strong shocks, if we substitute the shock condition (12) in (13), we get 

r*.,,,, .'ll.,,.,fx'i 


—+ (l/2)(/V-l) -2— dU^+ — 

y I ^ y y Pi) \y 


a„ (iV-l) + (A7Vr rf p„ 

~ 0- dr ^ NGmjXN)''^ dr _ 

p„ Un - \)+(XN)''^] r (A/- 1)+(XyV)*'* r’ ~ 


Substituting the values of da,, / ao,dpo/ Pq and m we find that 

—+5 —+ Cr''"' = 0 
dr r 


Where B = Afi.c = AC.A = —+ (i/2)(N - i)( 

Y N 
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B = 


2(N-l)y 


(N - \) + (XN) 


\n 


tv 


(x/y) 


2*,y(I - tv) pyv^///(l - w) 

1 

1 

n‘ 

p. p. p. 

\n - I)+ (XN)''^' 

r 


NGjXNyUnp^ 
[{N-\) + {XN)''^](3-w) 

On integration (30) disclose that 


U^ = kr-'’ 


B + { 2 -w) 


(31) 


Where K is constant of integration. 
The shock strength can be written as 


U 


\^(i / 


k'/ 


H>-(fl + 2) 


kUB + 2 -w) 


(32) 


Finally, the expression for the pressure, the density and the particle velocity 
just behind the shock for above cases can be expressed as 


Weak shock with weak magnetic fleld: 

2"'"’+T*|*r[p^2+3H')+2H'-4]/8 

P ~ KiK 


(33) 


l^tr-“+p,*f[p’(2+3w)+2>t'-20]/8 

M = kk.r'^ •' 


I* 


Weak shock with strong magnetic fleld: 

[l/pM2 + 3H')+6w-l2]/8 

p = + yA:, kr 

— f(l/p’)(2-.lB)+2B-2ol/* 

p = p + p, kr'- ‘ 

< ' 

Strong shock with srong magnetic field : 


(34) 
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P = 


!h. 

*2 


XiN) 


7-(h+«) 

kr 


(B + 2-w) 


+ y'(A/)*/“'■*’ 


p = N){kr-^ 


C 

(B + 2-w) 



(35) 


RESULTS AND DISCUSSION 

The expression (23) and (27), in general, represent the propagation of weak 
diverging spherical shock through a self-gravitating gas in the presence of weak and 
strong magnetic field, respectively. 

For — = 1.25 at r =100, the integration constant k and k for weak shocks in 

the presence of weak and strong magnetic field are given in tables 1 and 2 for w = 2 
and 2.5. Similarily for strong shock with strong magnetic field, when U/ = 20 at 
r = 100, the calculated values of c/k^ and k / kj are given in table 3 forw = 1.75,2.0 
and 2.5. Now calculation can be done easily. The conclusion is that the effect of the 
magnetic field (p^) is significant on the flow pattern. 


TABLE-1 


TABLE-2 


Weak shock with weak magnetic field 


W 

K 

0.01 


0.3979135 

0.10 

2.00 

0.2628988 

0.25 


0.1317615 

0.01 


0.0221839 

0.10 

2.50 

0.0135609 

0.25 


0.0059708 


Weak shock with weak magnetic field 


P^ 

W k 

1.75 

-5.2862155x10' 

2.00 

2.00 -0.0154822 

2.50 

-0.0442566 

1.75 

- 0.057393 

2.00 

2.50 -0.1787856 

2.50 

-0.5571584 
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TABLE-3 

Strong shock with strong magnetic field N=5, = 1.4 



W 

k 

C 

B 

2.00 

- 

792.39106 

-58.633031 

- 


1.75 



0.0465122 

2.50 


664.06213 

-68.345757 


2.00 


437.88152 

-53.914952 



2.00 



-0.1326276 

2.50 


459.00189 

-59.0754251 


2.00 


15.666711 

-105.41187 



2.50 



-0.2676282 

2.50 


18.093003 

-169.28974 
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DERIVATION OF A GENERAL CLASS OF 
BILATERAL FUNCTIONS FOR LAGUERRE 
POLYNOMIALS OF TWO VARIABLES 
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ABSTRACT 

A new class of bilateral generating functions for the Laguerre polynomials of two variables 
have been obtained. 

Mathematics Subject Classification (2000): 33C 45 


INTRODUCTION 

Applying Lie-group theoretic of W. Miller, we obtain the following operator for 
by giving suitable interpretation to n\ 




such that 


(a- , >’)/"] = -i- — L where ^ = t 

The extended form of transformation group generated by R is 




X . t . 

- -».v, - 

1 - wt 1 - wt 


where/(A, y, t) is an arbitrary differentiable function. 


( 1 ) 


( 2 ) 


( 3 ) 
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The object of the present is to derive a general class of bilateral generating function by 
employing group-theoretic method. Actually our result can be put in the form of a 
theorem as follows: 


THEOREM 

If there exists a unilateral generating functions of the form 


n = 0 

then the following of bilateral generating functions will hold: 


( 1-0 




X 




tz. 


(4) 



n*0 


a,, (pO'(P + «)! 

(p-«)!(n!)'(P + p)! 


(«.|3) 


(jc,y)r'’r" 


(5) 


Importance of the result of Eq. 5 is that whenever one knows a generating functions of 
the type given in Eq. 4 for a particular value of a then the corresponding bilateral 
generating functions can at once be written down from Eq. 5. Thus, one can derive a 
large number of bilateral generating functions by setting different values to a^. 


Let 


DERIVATION OF GENERATING FUNCTION 


G{x,y,t)= X 


/i = 0 


Now replacing t by tz, we have 


a 


(a.p) 


(x,y)r*z" 


( 6 ) 
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We operate both sides by exp (w^) 


exp (wi?) G (jc, y, tz) = exp(w^)^a,i'j“’^^(-«.)'Kz" 


Left hand side of Eq. 7 becomes 


(A) 


1-wt 


On the other hand, right hand side of Eq. 7 reduces to 


yy viz"'/" 


Equating (A) and (B) and substituting w = 1, we obtain Eq. 5. 


PARTICULAR CASE 


Setting^=0 and p = 0, we have the following bilateral generating function for 


m- 


0{x.l)-'^a,L“{xy 


then the following bilateral generating function will exist 


(1 - £ X “.f 

1 * A t J \ ^ / 


This was derived by .4/ salam. 
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ANALYSIS OF TRAFFIC NOISE PARAMETERS IN 
RELATION WITH PHYSIOLOGICAL EFFECTS 

V. K. KATIYAR*, K. S. BASAVARAJAPPA AND NILAM* 

(Received 12.12.2001) 

ABSTRACT 

Disturbance caused by noise is one of the most important health consequence. 
Over a number of years, investigations have shown that noise affects different activities 
and becomes the causative factor of environmental pollution. In view of the rapid 
infrastructural developments it is essential to study the traffic noise with respect to 
various causative factors (noise 65>dB leads to the annoyance i.e. above the 
permissible range). Keeping in view of the alarming increase of noise due to heavy 
traffic, engine sounds and blowing the horns, it is necessary to approximate the noise 
parameters. Present study takes into account the noise measurements on NH-58 at 
different noisy locations. Data is recorded using Noise Analyzer B&K2260 with 
software BZ7815. 

Keywords : Noise, traffic, sleep disturbances. 

INTRODUCTION 

Sounds are pressure waves travelling through air (or through fluids and solids). 
The waves, which are produced from the source point, will result in mechanical 
disturbances in material medium. Pressure variations propagate away from place of 
origin by transfer of energy from molecules are well associated with the pollution of 
dust particles which include chemicals (like carbon monoxide, nitrogen, hydrogen 
etc.). Sound pressure waves produced due to diesel generator vehicles moving on the 
road release the smoke and the sound (the engine sound and the horn sound) are the 
added parameters to the environmental pollution. Ultimately, the molecules, which 
carry quantum of these unwanted waves, cause the hazard on the human beings health. 
We notice the damage of hearing due to periodic compactions and rarefaction in the 
eardrum on tympanic membrane, the reduce of work efficiency the arrest of respiration 
due to the large quantity of carbon and dust particles enter the lungs and spech 

* Department of Mathematics, I.I.T. Roorkee, Roorkee, India; H-mail: vktmafma@iitr.emet.in 
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interference within the communication and other related problems. The regulation 
prohibiting the sounding of motor horn in a built-up area on any road is specified for 
the fixed time duration. 

Ahrlin and Ohrstrom [ 1 ] made a review of different kinds of noise exposure 
effects in relation to medical effects. Talbott et al. [6] focused on the study of 
occupational noise exposure for the increase of blood pressure levels. Van Dijk [7] 
reviewed the non-auditory effects of noise in industry comparing to stress factors. 
Griefahn [3] investigated the cardiac responses due to military noise at the time of 
sleep. Makarewicz [4] gave the model for A-weighted sound exposure level in a 
built-up area by investigating excess attenuation. Yoshida [9] experimentally 
demonstrated that noise at low or moderate levels could affect performance. Watts 
[8] considered the study of traffic running on concrete and bituminous road surfacing. 
Omiya et al. [5] examined the noise dose relationship and ratings of environmental 
noise by and with the help of stored data and concluded that noise environments 
irrespective of time of day and land use. Ando et al. [2] formulated for subjective 
attributes for noise fields based on the model of the auditory brain system, the sensitivity 
of human ear to sound source formed by physical system being compared with response 
of the cochlea. 

In this paper a study to approximate the noise parameters due to traffic noise 
measurements in relation to physiological effects when human beings are exposed to 
traffic noise for a longer period. 

FORMULATION 

Sound pressure waves have been propagated as plane waves with the quantity 
velocity, time, displacements, frequencies and the selected noise parameters. The 
noise can be assume to the unidirectional motion from the source point to the receiving 
point, the level of the noise varies in terms of various noise producing vehicles such as 
: noise from the engine, noise from the speed and noise by blowing the horns. Taking 
X as the displacement, t as the time, u as the velocity and c as the reference velocity, 

dx" 


(1) 
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U„^c^V^u(x,t) 


( 2 ) 


where 


2=^+_?L+-£L 

dx^ dy^ dz^ 


(3) 


Equation (3) is employed in equation (2) for one-dimensional system. The ranges of 
'x' and 7* are considered at three different locations and timings for monitoring the 
noise measurements. Measurements are rq)eated in different time steps also. Distances 
for monitoring are taken at 5,10 and IS meters from the center (noise source point) 
of the road. The effective value of the noise monitoring is the best measure for 
continuous sound at short and long duration. This assessment is characterized by a 
parameter called crest factor which is the peak sound pressure level to the RMS 
value of the sound. It is ^proximated as the logarithmic ratio of peak levels and RMS 
values with reference to the sound wave. 


C,=K,[K,,-K,,] 


(4) 


where, K, = 20, AT,, = logp^,, K,, = logp^ 

ANALYSIS 

Solving the equation (1) by the closed form technique (tentative solution with 
specification of variables) as 

u{x,t) = X{x)T(t) (5) 

where, A'is a function of x alone, 
r is a function of t alone. 

Introducing (5) in (1) we get after discritization 


T"{>) 

m XU) 




( 6 ) 


The value of k is for three different variations k>0, k = 0, k <0. 

From equation (6) taking separate differential equations in t and x respectively and 
solving them we obtain the solutions with successive arbitrary constants. The solutions 
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for /c > 0 and A = 0 do not represent the physical nature of the problem. Since the 
terms with k<0 involve the way parameters. Therefore, the solution of (1) for ^ < 0 is 


u{x,t) = ^Acos^fia +Bsinyfia^ 


Ccod 


\fk 


\ ^ J 


Jc + Dsid 


Jk 


( 7 ) 


The constants A, B,C,D have been evaluated by prescribing the boundary conditions 
with variations of displacements x from 2 to 15 meters and time t from 1 to 30 hours. 

The term yfk is replaced by 2nf per second (the frequency function which takes the 
noise measurements for showing the values of L^. Values of u (x,t) have been calculated 
for the noise measurements taken at three locations (shown in the table 1). Noise 
measurements have been recorded at three locations for different types of vehicle 
sound and the noise values have been employed in equation (7). 

In terms of energy emission levels from the source point to the receiving point it 
is necessaiy to reduce the noise levels in terms of various adjustments such as allowing 
the vehicles in the lanes adjusting the reaching distance alternately, the barriers with 
different heights. 

Traffic flow adjustment (AT) 


AT = 101og(I^7iDo/5,T) (8) 

It is necessary to analysis the adjustment of traffic flow in order to reduce the noise 
effect in term of distances for vehicle into the traffic stream. This is accomplished 
through the application of the traffic flow adjustment. Three of the five subparameters 
in this adjustment are constants, with n representing an integration constant for the 
assumed infinite roadway segment (acoustic angle of influence is n, or 180 degree). 
D,, represents the reference distance of 10m, and T, the time period, is 1 hr. 

10Iog2I^. = 10(0.3) +lOlogF, (9) 

Conversely, if speed was doubled and everything else held constant, the result would 
be a 3-dB decrease in Leq. 


10 log 


_1_ 

25; 


lOlog^ + lOlog 


]_ 

5 , 
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= -3 + 101ogJ- (JO) 

The decrease in AT is that, with volume held constant doubling the vehicular 
speed actually the amount of time that each vehicle acoustically influences a given 
receiver of sound. 

Distance adjustment (AD) 


AZ) = 101og(Do/D) 

The distance adjustment or spatial decay rate can be established to be 3 dB per 
distance doubling, ifZ) is doubled. 


101og^ = lOlog^ + lOlog 


_ 1 _ 

D 


= -3 + 101og^ 


( 11 ) 


Barrier adjustment (AB) 

The placement of an obstruction between the highway and the observer will 
shield the observer. 

In most instances where the shielding adjustment is used, solid barriers or walls 
are involved. Barrier attenuation analytically defines the amount of acoustic energy 
loss encountered when sound rays are required to travel over and around a barrier. 

The Fresnel number may be calculated as 

iV = 2| (12) 

where, 5 is the path length difference, X wavelength. 

Thus, a barrier with 8 = 0.9144wj would have A^=273. 

RESULTS AND DISCUSSION 

The elevated noise measurements measured at noisy locations describe the 
traffic noise effects on various geometrical dimensions keeping in view that the 
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allowable range is 65 dB. The Frequency range from 800 to 4000 hertz is 
employed in the analysis. It is clear that the road traffic noise data measured 
disagree with the allowable range. These elevated values of noise ensure the 
sleep disturbance, activity interference for the human beings living close to the 
noise source. Also the prolonged exposure with these elevated noise cause the 
hearing loss and damage on the hearing system which becomes the significant 
physiological effect. The angular velocity of the sound wave varies from the 
allowable (300 m/sec) to the elevated level (i.e.>400 m/sec) justifies the 
corresponding increase in the normal velocity. From equation (7), the increase 
of velocity by 11.5184 m/sec induce the temporary threshold shift at frequency 
3500 Hz and further increase of velocity by 13.90696 m/sec will become the 
cause for the effective threshold shift indicating physiological destruction within 
in the hearing mechanism. The higher range of crest factor C^46.1-52.26 dB at 
5 to 15 m show the best measures (equation 4) to analysis the short and long 
sound waves. As the distance from the source increase there is a decrease in 
crest factor, this shows that as the distance from the source increases, effect 
ofnoise diminishes. It can be concluded that some noise control techniques may 
be employed at noise source transmission path and at the receiver. At source, 
reducing the exciting forces and changes in operating procedure can control noise. 
At the transmission path, it can be controlled by the construction of the barrier 
(equation 8-12) and absorption devices. The receiver can use the high quality 
protective equipment in order to avoid consequences of the noise. 


Table 1 : Hourly Variation of Leq. 


Location 1 

Location 2 

Location 3 

Time 

Leq 

Time 

Leq 

Time 

Leq 

10-11 

80.2 

. 10-11 

82.3 

10-11 

77.3 

11-12 

83.1 

11-12 

81.6 

11-12 

77.7 

12-1 

82.2 

12-1 

81.5 

12-1 

79.1 

1-2 

82.3 

1-2 

82.6 

1-2 

77.5 

2-3 

82.2 

2-3 

81.8 

2-3 

78.6 

3-4 

78.7 

3-4 

81.9 

3-4 


4-5 

78.9 

4-5 

81.7 

4-5 

77.9 

5-6 

82.1 

5-6 

82.4 

5-6 

79.2 
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Table 2 : Calculation of Crest Factor (Cp. 


Displacement 

LA^^ (peak) 

c, 

Cj (Average) 


105.0 

52.474 


5 meters 

105.1 

52.472 

52.56 


105.1 

52.708 



97.5 

45.656 



102.2 

47.467 


10 meters 

102.2 

49.218 

49.69 


103.6 

50.200 



103.6 

49.957 



102.3 

50.445 


15 meters 

97.9 

53.452 

46.10 


97.9 

44.461 



98.1 

46.510 



Table 3 : Classified Trafflc Volume at Various Locations. 


Location 1 


Time 

C/J/V* 

Scooter 

M/C** 

Bus 

Truck 

Tempo 

Tractor 

lOtoll 

205 

193 

107 

81 

46 

44 

25 

11 to 12 

187 

203 

97 

87 

37 

38 

22 

12 to 1 

193 

249 

84 

80 

32 

34 

36 

1 to 2 

275 

369 

132 

86 

42 

11 

19 

2 to 3 

175 

187 

104 

80 

41 

20 

19 

3 to 4 

173 

97 

59 

79 

42 

17 


4 to 5 

177 

130 

40 

61 

41 

26 

19 

5 to 6 

181 

175 

HI 

85 

47 

31 



Location 2 


Time 

C/J/V 

Scooter 


Bus 

Truck 

Tempo 

Tractor 

lOtoll 

211 

163 

120 

154 

79 

112 

67 

11 to 12 

191 

154 

99 

157 

88 

122 

57 

12 to 1 

201 

140 

110 


61 

75 

26 

1 to 2 

214 

146 

52 

126 

63 

61 

37 

2 to 3 

140 

201 

88 

128 

72 

68 

40 

3 to 4 

177 

212 

113 

105 

51 

90 

39 

4 to 5 

270 

336 

105 

120 

102 

78 

27 

5 to 6 

245 

259 

117 

111 

92 

81 

21 
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Location 3 


Time 

C/J/V* 

Scooter 

■SSSil 

Bus 


Tempo 

Tractor 

10 to 11 

190 

165 

74 

41 

65 

87 

40 

11 to 12 

168 

196 

80 

35 

72 

88 

44 

12 to 1 

156 

100 

45 

56 

21 

88 

32 

1 to 2 

122 

58 

36 

42 

44 

51 

31 

2 to 3 

154 

71 

38 

64 

58 

66 

27 

3 to 4 

174 

125 

66 

36 

78 

81 

35 

4 to 5 

182 

142 

84 

76 

70 

81 

39 

5 to 6 

176 

121 

90 

56 

65 

82 

25 


*C/J/V-Car/Jeep/Van 
** M/C - Motor Cycle 


Table 4 : Individual Vehicle Noise Levels. 


Vehicle 

Type 

Car 

Jeep 

Van 

Scooter 

Motor 

Cycle 

Tempo 

Bus 

Truck 

Tractor 

Noise 

Level 

(dB) 

76.14 

78.01 

77.80 

75.10 

76.55 

84.14 

86.1 

87.30 

89.9 


Velocity Vs Distance 



DIabincu (m) 


Fig. 1 : Velocity Vs Distance 
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Fig. 2 : Spectrum Showing Leq 
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ABSTRACT 

In this note, we give three decompositions of continuity which are not widely known. For example 
every text-book in calculus contains following: (a) A continuous function preserves compact sets, 
(b) A continuous function preserves connected sets (the intermediate value theorem). That (a) and 
(b) together give a nice space is not found in text-books. 

Mathematics Subject Classifications (2000): 54 C 05,54C08. 

Keywords and Phrases: Continuity, connected function, compactness preserving function, 
decon^ositions. 

1. INTRODUCTION 

In the mathematical literature there are concepts and results which are simple, yet 
profound and beautiful. Some of them can be understood by undergraduates but are 
not widely known. In this note it is proposed to present a few gems related to continuity. 
The results are from elementary analysis and have successive generalizations to metric 
spaces, topological spaces, and topological vector spaces. In order to make this note 
tangiable, we give proofs for functions on IR to IR. Then in the remarks we suggest 
how the generalizations are conceived. We provide a reasonably complete bibliography 
which will be useful to those who wish to pursue the subject further. 

The main theme in the present note is the decomposition of continuity into two weaker 
conditions. The classical case that continuity is equivalent to upper and lower semi¬ 
continuity is widely known. Here we present three other such decompositions. The 
concept of continuity being one of the most important concepts in mathematics, any 
such decomposition will naturally shed more light on it. 

2. COMPACT AND CONNECTED 
Suppose /://?->//? is continuous. The following results are well known: 

(a) / takes compact sets into compact sets, and 

(b) / takes cormected sets into connected sets (Intermediate Value Theorem). 

But it is not widely known that (a) and (b) together imply that / is continuous! (See 
[ 10 , 171 ). 

* Emeritus Professor of Mathematics, Lakehead University, 96 Dcwson Street, Toronto, Ontario. M6H, 1 H3, Canada 
Website: wvvw.naimpally.com E-mail: sudha@accglobal.net 




40 


WHEN IS A FUNCTION CONTINUOUS 7 


(2.1) THEOREM. If f :IR-^ IR satisfies (a) and (b) then it is continuous. 

Proof : It is sufficient to show that / is continuous at 0. Suppose it is not. Let 
U„ = (-« ' ,«■'). Then there is a closed interval with/(;c) as midpoint such that for 
each n eIN,f{U„)iBQ. Thus there is an x^ such that /(x,) ^B^ Hence there 
is a closed interval 5, with midpoint /(jc,) and disjoint from B^. Since/preserves 
connected sets, /(C/^)czand so there is an such that 

u 5,. Then there is a closed interval B^ with midpoint /{x ^) and disjoint 

from u 5,. Inductively we have e U„ with /(jc„ )« .0 < ^ n -1}, fi* a 

closed interval with midpoint f{x^) and {B/^ikelN] pairwise disjoint. Clearly 

;c^ -> 0. So the set E = {0}u{jrj:A elN} is compact but f{E) is not compact, a 
contradiction of (a). 

(2.2) REMARK. Suppose X and Y are metric spaces and f:X ->Y satisfies (a) 

and (b). Then, noting that in the above proof each U„ is connected, we find that /is 
continuous if A" is locally connected ([ 10]). This result was extended in [WH] to 
topological spaces where A'is locally connected and first countable and Y regular. 
The paper contains counter examples to illustrate the need for these conditions. After 
about a decade, it was shown in [7] that White's result is true if Tis merely Hausdorff 
rather than regular. 

3. CONNECTED AND CLOSED FIBERS 

If f :lR-> IR is continuous, then for each y elR, the fiber /"‘(y) closed and, of 
course,/preserves connected sets. The following converse was proved in [16]: 

(3.1) THEOREM. Suppose /://?->//? takes connected sets into connected sets 
and for each y elR, the fiber /"' (y) is closed. Then/is continuous. 

Proof : Suppose x e/R and /(x) ^{c,d). Then f~'\c,d] is a closed subset not 
containing X So there is an open interval I containingx and disjoint from /"' {c,d}. 
Since/(1) is connected and contains f{x),f{I) c {c,d ). Thus/is continuous. 

(3.2) REMARKS. To understand how this result can be generalized, we note the 
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following: 

(a) Open intervals are connected and form a base for the topology of IR. 

(b) The set {c,d] is the boundary of the open interval (c.d). 

Rowe's result (3.1) was successively generalized to metric spaces [10], to Hausdorff 
topological spaces [ 12] and finally to a topological result with qjplications to Functional 
Analysis [8]. We give below Hrycay's result whose proof is as easy as that of (3.1): 

(3.3) THEOREM. Suppose is a topological space with a base p of connected 
open sets and K is a regular topologicl space with a base p of open sets. Suppose 
/; A' y is a function satisfying; 

(a) for each R e P, f{clB) is connected, (i.e./is connected- p), and 

(b) for each 5'€P' , / "' [3 fl '] is closed. (55' is the boundary of 5'.) 

Then/is continuous. (See [8].) 

The result (3.3) includes, as special cases, most of the previous results of this type 
and moreover, has applications to Functional Analysis. In fact, Hrycay arrived at this 
result in an attempt to solve a problem in Functional Analysis. In 1966,1 was teaching 
a course in Functional Analysis and also conducting a seminar in connected and other 
non-continuous functions. In the course we studied the result: 

(3.4) A Banach space A"is finite dimensional if any and if every linear functional on it 
is continuous. 

It is reasonable to conjecture that if the Banach space A^is infinite dimensional, then 
any linear functional on A, although not continuous, ought to satisfy some weaker 
property implied by continuity, such as preservation of compactness or connectedness. 
So I made the following conjecture: 

(3.5) CONJECTURE. Every linear functional on a Banach space preserves 
connected sets. 

Hrycay who was looking for a topic for his Thesis, took up this problem and showed 
the following : 

(3.6) For every topological vector space (TVS) A, there is a base p of connected 
open sets such that every linear functional on Ais connected- p (see (3.3) (a).) 

From (3.3) and (3.6) he deduced the well known result: a linear functional on a TVS 
is continuous if and only if the fibers are closed. 

(3.7) A linear functional on a TVS is continuous if and only if it takes connected sets 
into connected sets. 
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Thus the conjecture (3.5) is false. Hiycays woik showed die significance of the generalization 
from connected functions to functions that are connected- p. For details see [8]. 

4. QUASI AND NEAR CONTINUITY 

Two weaker forms of continuity have been in the literature for over half a century and 
are widely known to workers in Topology and Functional Analysis. Before describing 
them, let us recall that 

(4.1) a function f:X -^Y is continuous at x eX if for each neighbourhood V of 
fix), there is a neighbourhood U of;c such that U c. 

By changing the condition "U c /‘'(F)" in different ways, we get weaker forms of 
continuity. 

(4.2) DEFINITION. A function f:X -¥Y is nearly continuous at xe^X if for each 
neighbourhood V of fx), there is a neighbourhood Voix such that U c 

We say that /is nearly continuous if it is nearly continuous at each xeX. 

Obviously a continuous function is nearly continuous and the following well known 
example shows that the converse is not true: g:IR-> IR given by 

g(A:) = 0 when x is rational, and 
= 1 when X is irrational. 

(4.3) REMARKS. Blumberg [2] discovered nearly continuous functions which he 
called "densely approaching", certainly a more appropriate term. Some authors have 
called such functions "almost continuous" which caused some confusi(m since some 
other weaker forms of continuous functions are also known under this term. The 
importance of near continuity stems fi'om the result: 

(4.4) Every linear function from one Banach space to another is nearly continuous. 

Thus near continuity and its dual, near openness, form intermediate steps in the well 
known closed graph and open-mapping theorems in Functional Analysis. Moreover, 
Blumberg proved a spectacular result: 

(4.5) Every function /;//?-►//? is nearly continuous on a dense subset of IR. 
(i.e./ has the Blumberg property). 

There is a vast literature on the Blumberg property (see [ 13]). 
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(4.6) DEFINITION. A function f:X is said to be quasi-continuous at xeX 

if for each neighbourhood Uofx and each neighbourhood V of/(x), there is a non¬ 
empty open set fVaUn A quasi-continuous function is one that is 
quasi-continuous at each x eX. 

Clearly if fV contains x, then/ is continuous and so continuity implies quasi-continuity. 
That the converse is not true is given by the following example; h: IR -> IR given by: 

h{x) = sin(jc'') when x 0»and h{x) = 0 when x = 0. 

(4.7) REMARKS. The following is a well known example of a 
function s:IRx IR-> IR that is separately continuous but not jointly continuous: 

5(x,y) = xy{^x^ when(x,>') * (0,0) 

= 0 when (x,>') = (0,0) 

It was first observed by Volterra, who was quoted by Baire [ 1 ], that 

(a) every sqjarately continuous function f : JR x IRIR is quasi-continuous. Moreover, 

(b) every such quasi-continuous function is continuos on a residual subset of IRxIR. 

From (a) and (b) follows a classical result due to Baire that separate continuity of 
f.IRx IR-y IR impliesthejointcontinuityof/bnaresidualsubsetof//?x //?. Thus quasi¬ 
continuity serves as an important intermediate step in the proof of Baire's result. There is a 
vast literature on this topic of separate versus joint continuity (see 11,14,15). Quasi¬ 
continuity was systematically sudied by Kempisty [9] and has been rediscovered [3 ]. 

Although these two important generalizations of continuity viz. near and quasi continuity 
have been studied extensively for a long time, it was only relatively recently recognized 
that they give a nice decomposition of continuity (see [4,13]). 

(4.8) THEOREM. If f:IR->IR is nearly continuous and quasi-continuous, then it is 
ccHitinuous. 

Proof; Suppose/is not continuous at x e IR. Then there is a closed interval B, with 
midpoint f(x), such that in every open interval I containingx, there is ay with f{y) ^ B . 
Since/is quasi-continuous aty, there is an open set IF c / such that /( W)ciB^. Then 
Ictclf'\B) which contradicts the fact that/is nearly continuous. 

(4.9) REMARKS. The above result is a global one and does not hold locally. For 
example the function j:IR-* IR given by y^x) = sin(x’') when x is rational, and 
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j{x) = 0 when x is irrational is nearly and quasi continuous at 0 but is discontinuous 
every where. In case of topological spaces, a function that is nearly and quasi continuous 

is 0 - continuous. A function f:X -^Y is 0 - continuous iff for each x eX, and 
each neighbourhood Voff (x), there is a neighbourhood Uofx such that f{U)czclV. 
Clearly if Kis regular, then 0 - continuity is equivalent to continuity. We have discussed 
here three decompositions of continuity. There are others; See e.g. [5,6]. 
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ABSTRACT 

In this paper values of Pie determined by various mathematicians from ancient to modem 
period are given. Also we present a number of series representation for Pie. 

Keyrvords : Tritah; Circumference; Diameter; Circle. 
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Many mathematicians studied the values of Pie from very beginning. Some of 
them are Aryabhata I (b. 476. A.D.) Brahamagupta (b. 598 A.D.), MahdviracSrya 
(b.629 A.D.) Bhaskaracdrya, Ganes'a Daivajna (b.l507 A.D.) [1-16, 20, 23, 24, 
27]. 

The oldest value of Pie is found in Rgveda (2500 B.C.) (for instance, [8, 24]) 
According to Rgveda.- 

Bhinad Valasya Paridhin iva tritah (Rgveda 1-52, 5; [8], p.l87) 

This means 

The ratio of circumference (paridhi) and diameter (vyds) is 3:1 and it is called tritah. 
That is, 

— = — = tritah 
d 1 

Where c and d stand for circumference and diameter of a circle respectively. 
Baudhdyana (1000 B.C.) suggested that value of Pie is equal to 3 [39]. 
Archimedes (b.287 B.C.) had shown the Pie to be the ratio between 

223 

circumference of a circle and its diameter, and found that it lies between and 
220 

[18-20,23]. Archimedes was the first who introduced the notation ji [40]. 
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After Archimedes, Liz Hui (Chinese mathematician(c.263 A.D.)) calculated 
the gross value of Pie = 157/50 = 3.14 and near value 3927/1250 = 3.1416 [23]. Tsu‘ 
Chung-Lhih (b. 470 A.D.) in Tru Su gave the value of Pie equal to 355/113 [23,24]. 
For the details see; Ancient Chinese Technology and Science, pp. 52-123. 

For detailed sources on the simplest approximation 3, refer to [26]. These 
include sources such as Jaina, Buddhistic, Christian, Chinese, apart from Hindu 
sources. Mishra and Singh [26] also present a good account of the 


22 355 600 16 n 6 \ 

values—>V1",-,-,— and — 

7 113 191 5 I 9 j 


for Pie in almost all culture areas of the 


world. Indian knowledge for Pie is further enriched by [25] where authors establish 

355 62832 

Virasena's knowledge for 777 >n> 7 --— through first degree indeterminate 


113 


20000 


analysis concept in the equation c = 3d + of Virasena (a ninth century 


113 


A.D. Jaina mathematician). 


Aryabhata 


/(b. 476 A.D.) was the first Indian mathematician who gave the value 


of Pie in his Aryabhatiam (c. 499 A.D.) correct to four decimal places [6,10,15, 
21,23,24,28, 34]. He says the following verses: 


I ([24], p.28) 

This means 

100 plus 4 multiplied by 8 and added to 62000, this is approximate measure of 
circumference of circle whose diameter is 20000 . 

c (100+ 4).8 +62000 62832 
d 20000 20000 


After Aryabhata I, Brahmagupta md Mahaviracarya considering many figures 
concluded that circumference divided by diameter is 3 or square root of 10 [15,21, 
23,34]. In other words, for concluding the value of Pie they could not be as accurate 
as their predecessors. 

Bhaskaracarya //(b. 1114 A.D.) gave the formula in his LilavSti for Pie [1,3, 6 , 
9,11,22,23,27] in the following stanza. 
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czn^ ^ ^CT: I 

glf4?lf^ ir f^SSI M: We['«|it|gK Tfnr^i | 

([27],p.l43) 

This means 

Diameter of a circle multiplied by 3927 and divided by 1250 gives the near 
circumerence, alternatively, the diameter multiplied by 22 and divided by 7 gives 

3927 22 

the gross value. Near value of Pie is equal to and gross value —. 


After Bhaskaracarya II majority of mathematical work was in progress in 
southern part of India, and Kerala's mathematicians determined the series for Pie. 

Patuman Som ayaji in his Karanpaddhati (c.l438 A.D.) gave the series [28, 
35] 


In Tahtra Sahgraha (c. 1502 A.D.), NUkanfha had been anticipated this series 
before Gregory. 

After two hundred fifty years this very series was also enunciated by European 
James Gregory (1638-1675 A.D.) in 1671 [23,24,28]. 

* -I 

tan x = x -+- 00 

3 5 

Ifwetake;c= 1, 


or 


4 

c = 4rf- 


i 1_ 

3*^5 
4rf 4rf 
3 5 


00 


00 


From this series Pie can be found for ten decimal places, i.e. 

71 = 3.1415926536 

The Kerala School provided many approximations to Pie. One such value of 
Pie is given in Kafapayddi notation [28, p. 177] is that the circumference is 
represented by 



and diameter by 








48 


A HISTORY OF DETERMINING VALUES OF n FROM ANCIENT TO MODERN PERIOD 


Then the good approximation of Pie = 


31415926536 


10,000,000,000 

= 3.1415926536 which is correct to ten decimal places 

Madhava (1340-1425 A.D.) gave rational ^proximation for Pie by considering 
finite number of terms [11,12,15-17,28]. 


n = 


'll ^ I . K 

1 1-i-1-; 

3 5 2K-1 4K^ + l 


Alternatively, above series for K = (/i+l)/2 or n = 2K-1 may be expressed as 
n = 4 ! 


1.14- 

3 5 n (/i + l)'+l 


Last but one term is negative or positive according as is even or odd. 

If we take « = 1, Pie comes out to be 3. 141564 which is correct to four decimal 
places, 

Other rational approximation given hy Madhava is 


^ - 1 ^ ^ 


±-T 

n 


n + \ 


+ 1 


- 

fn + V 

U1 

t/I + ll 

4 

1 2 , 


2 


He also gave infinite series 
7t = Vn[ 


i-L.. ‘ 


1 


-00 


3.3 5.3^ 7.3^ 

From the above series Madhava gave a rational approximation to Pie as [28] 
2827433388253 


7t = 


9 " 


= 3.14159265359 
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Which is correct to eleven decimal places. 

Nilkantha Somayaji (b. 1444 A.D.) in his Tahtra Sahgraha gave some series 
expansions for Pie [28,35,37]: 

71 1 1 1 

i”2^-l'^6'-l'^10'-l 
1 1 1 

=- 1 -- 

1.3 5.7 9.11 

4-7I_ 1 1 1 

8 ~4'-l’^8'-l'^12'-l'^ 

£-3^J _J_1_ 

4 ” 2.3.4 4.5.6 6.7.8 8.9.10 

71-3 1 _1_ 1 

6 (2.2^ -1)^ - 2^ ^ (2.4^ -1) - 4^ ( 2 . 6 ' -1/ - 6^ 

_ 7^__1 _ 1 _ __1 _ 

16” 1’+4.1 3’+4.3 "^5*+ 4.5 

German mathmetician Vieta (b. 1540 A.D.) showed that Pie can be found from the 
following formula[20]: 

2 _ IT ipr IT h j_ nr_ 

n~n\2'^2nu'*'2U2 

English mathematician John Wallis (1616-1703 A.D.) also showed a number pattern, 
which was entirely new way of closing in on the value of Pie [20] 


71 2.2.4.4.6.6.8.8- 

2 1.3.35.5.7.7.9.9 
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Lord Brouncher searched for a formula for Pie [20]. 


n - 


2 + 


25 


2 + 


49 


81 


2 + 


81 


121+WC. 

Sankara Verma (c. 1819A.D.) in Sadrathnmdld [28] gave Pie as 


314159263358979324 ^ 3 i4i5926558979324 
10 ” 

which is correct to seventeen decimal places. 

Lindemann (b.l852 A.D. in 1882 is credited for proving the transcendence of Pie 
[13,20,24], 

SvamiBhdrti Krsna Tirtha (1884-1960 A.D.) the Sankardcarya of Govardhan Math, 
gave Pie in Katapayddi notation [2, 28]. To know more about Katapayadi system 
[27, pp. 175-176] 

11 ([2],) p. 320) 

This gives the Pie for thirty one decimal places, that is, 

71 = 3.1415926535897932384626433832792- 

Ramanujan (b, 1887) also obtained some approximations for Pie [28-31] 

Some finite expansions for Pie are: 


—V7 = 3.14180. 
16 


88 


7-3V2, 


= 3.14159274. 
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355 0003 ^ 

113 V 5533 J 


3.1415926535897943. 


He also gave infinite expansions for — [28,30] 

7C 


n 64[2) 64V2.4J 


( 1 ) 


1 ^ 1103 27493 1 1.3 53883 1.3 1.3.57 _ 

27tV2 “ 99' 99" 2 3'"^ 99'" 2.4 4'8' 

In series (2), 

1103 

-^ = 0.11253953951, 
while ^ = 0.11253953951, 

which is correct to eight decimal places even by taking only I term. 

This very series (2) was used in 1986 to compute Pie for 17 million places 
{Ganita Bhdrati 1994, 84-85). 

In 1998, Jasumasa Kanad (Jj^anese) found the Pie for 20,1326000 decimal 
places [24]. 

In 1989 David Chudanovaski and Grary Gregory of Clombia University found 
the pie for 1011196691 decimal places [38]. Surprisingly, in 1995 Takahashi and 
Kanada obtained 3221225466,4294967286 and 6442450938 digits for Pie which 
were extended by them in 1997 to 51539600000 digits. 
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MAGNETIC BEHAVIOUR OF TERBIUM-IRON- 
TRANSITION METAL MIXED OXIDES 

U. P. SINGH*, A. N. THAKUR* AND V. SINGH* 

(Received 12-08-2002 and in revised form 30-11-2002) 

ABSTRACT 

The compounds of Terbium-iron-transition metal mixed oxides of the type WFeTO^ have been 
prepared by solid state reaction technique and characterised by XRD pattern. The molar magnetic 
susceptibility (Xm ) of the powdered san^le have been reported in the ten^rature range 300K to 
1100. All studied materials show a typical ferrimagnetic behaviour and xJ vs f plot can be expressed 
in terms of standard relations for.fenimagnetic materials. The slope of asymptotic line to the curve, 
yield an average magneton number which indicates that all materials are perfectly ionic. The transition 
metal ions are Fe^* and Fe^* in TbFefi^, F^* and Cr^ in TbFeCrO^, Fe^* and Mn^* in TbFeMnO^, and 
F^* and Co'* in TbFeCoO^,. 

Keywords: Magnetic susceptibility, Terbium-iron-transition metal mixed oxides. 

Physics Subject Classiflcation: 07.55 - W. 75.30 Cr 

INTRODUCTION 

Rare-earh and transition elements are characterised by their unfilled 4f and 3d orbitals 
respectively. Due to this reason, they yield variety of materials with intoesting and useful 
magnetic properties which have not only enriched our understanding but have found 
numerous application [1-3], We have prepared a series of rare-earth transition metal 
compounds with general formula R^JTO^ where stands for rare earth, T and T for 
transition metals. We have already reported the electircal transport studies of YFe TO^ 
and where T=Fe, Cr, Mn and Co [4-5] and magnetic behaviour of GdFeTO^ 

and YFeTO^ [6-7]. This paper reports results of our study on the magnetic behaviour of 
compounds with general formula TbFeTO^ with T=Fe, Cr, Mn and Co. So far, no study 
of this kind has been reported on these compounds. 
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MATERIAL PREPARATION AND CHARACTERIZATION 

The starting materials for the preparation of these compounds were Tbfi^, Fefiy, 
MnO^ and CoO which are 99.9% pure. The stoichiometric amount of these oxides 
were mixed and heated in a silica crucible for 50-hours at temperature of about 1400K. 
In this process, the mixture was subjected to one intermediate grinding and the final 
product was cooled down at slow rate. The prepared compounds under go the following 
solid state reaction. 


ITb.O, + 4Fe^O, 4TbFe^O, + 

Air 

imo, + 2Fe,Oy + 2C/-,0, > AThFeCr O, + Oj 

Air 

Tb,0, + Fe^O, + 2MnO. — 2TbFeMn 0, + 0, 


Tb,0, + Fe.O, + 2CoO ITbFe CoO, 

Air 

The weight loss corresponding to loss of oxygen of the right hand side of the 
above reaction was observed in all cases execpt in WFeMnO,.In this case, the 
observed loss was slightly less than expected. The details are described else where[8]. 

To confinn the complete formation of the prepared compounds XRD pattern 
were obtained for each material using CuKa radiation with X = 0.15418 nm. From 
XRD pattern the values of interplaner spacing £f^^,have been obtained using relation. 


_ 0.15418 


2 sin0 


( 1 ) 


From values, structure of the prepared compounds have been resolved using 

standard procedure[8]. All the peaks have been identified and assigned proper 
values. This confirms that prepared compounds are in single phase. All the compounds 
have been found to have orthorhombic unit cell with unit cell parameters a^, and 

as given in table 1. 
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Table 1. Structural parameters of the studied TbFeTO^ compounds 


Compounds Lattice Parameters (nm) Unit Cell Calculated density 




*0 

<^0 



TbFefi^ 

0.6246 

0.7366 

0.8836 

Orthoihombic 

5.34 

TbFeCrO, 

0.6340 

0.7482 

0.8850 

Orthorhombic 

5.14 

TbFeMnO. 

4 

0.6174 

0.7286 

0.8850 

Ortiiorhombic 

5.05 

TbFeCoO, 

4 

0.6280 

0.7410 

0.8915 

Orthorhombic 

5.20 


MEASUREMENT OF MAGNETIC SUSCEPTIBILITY 

Magnetic susceptibility measurement has been done using Faraday method. In 
this method one requires only 50-100 mg of the sample as well as the standard materials. 
The set up consists of a magnetic balance (of 10'^ gm accuracy with all control external, 
from Keroy, India), a tapered pole pieces (2" diameter from 4" base at angle of 30°) 
electromagnet and constant current power supply (both from Pollytronic, India). The 
pole gap is variable but has been fixed at a separation of 4.8 cm. At this pole gap, the 
magnetic field in constant //(J///cfe) region, which lies about 2.5 cm above the center 
of the pole is about 1.5x10* Am' at 2 A current in the electromagnet coil. A pyrex 
tube of 2 cm long and 0.5 cm dia has been used as a sample holder. Gd^ (IFO^)^ with 
molar magnetic susceptibility value 6.85 x 1 O’’ at 300 K has been used as standard 
substance. 

The standard material and the sample is hanged from the hook provided in the 
pan of the balance in constant H {dH/dz) region in the pyrex tube and weight is measured 
in both cases, with and without applied magnetic field. The molar magnetic susceptibility 
of the sample have been obtained from the relation [8,9]. 

= (Am / Amj )(mj / m)x^ 

Awi and A/«,are changes in the weights of the sample and standard substance. 
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m and are their masses and Xs the molar magnetic susceptibility of the standard 
substance. Being a relative method most of the error are automatically eliminated except 
the error in the measurement of mass m and weight change Am. The maximum probable 
error in these measurements has been about 2% at lower temperature (7’<500 K). But 
it increases with the increase of temperature, because hot air movement disturbs the 
sample holder in spite of closed one end of the furnace, and becomes as high as 5% 
around lOOOK [8,10]. 


RESULTS AND DISCUSSION 

The molar magnetic susceptibility (Xm) ofall the componds have been measured 
in both heating and cooling cycles. No hysteresis was observed in x„, and values were 
found to be almost same in both heating and cooling cycles. However a small weight loss 
was noticed in heating cycle probably due to presence of moisture. The results are presented 
in figure 1 as Xm vs T plot. It is seen from this figure that nature of all these plots are 
essentially similar. In general, 5 ^’' vs T plots are linear at higher side of temperature. 
However, there is systematic trends of experimental points towards temperature axis at 
lower side of temperature. The curves are thus similar to a standard ferrimagnetic matoial 
and systematic downward trend is due to the onset of short range magnetic interaction at 
lower side of temperature. We have tried to fit the experimental points to the standard 
equation of ferrimagnetism [ 11 ] given below by choosing the suitable parameters. 

_? 1 _ 

X. c„ c„(7--e) 

where c*, is the average value of Curie constant, 0„ is the asymptotic Curie temperature 
and and 0 are parametric temperature. The curves using eq. 2 are drawn by full line in 

respective vs T plots. It is seen that experimental point can be well fitted by eq. 2 over 

wide temperature range. The values of Cu, 0a,0* and 0 are given in table 2. The 
ferrimagnetic Curie (or Neel) temperature have been evaluated using condition 
7 ;, Xm' 0. This gives. 
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(7’.-ej(n-0) = 0?, (3) 

The real and positive value of T are meaningful and have been evaluated using above 
relation. These values are also given in table 2. 


Table 2. Magnetic parameters of studied compounds 


Compounds 

(K) 

0 

(K) 

e. 

(K) 

T C~ 

(K) (mor'm^k)x\0* 

TbFe,0, 

2 4 

-348 

236 

85 

252 

2.333 

TbFeCrO, 

-937 

243 

116 

255 

2.050 

TbFeMnO, 

-990 

165 

194 

198 

2.167 

TbFeCoO, 

4 

-540 

245 

134 

267 

2.250 


The compounds TbFeTO^ contains only three types of magnetic ions 
or F^* and or Thus at temperature much higher to T, the molar magnetic 
susceptibility of these compounds can be expressed by the relation. 

T'-e./r-e./r-e., 

where N is the Avogadro number, p p is Bohr magneton, is the permeability constant, 

k is Boltzman constant, P,, P^ and Pj are the magneton numbers of three types of magnetic 
ions respetively and 0ai,0„2 and 0„j are paramagnetic Curie temperatures which takes 
into account the effect of various interaction. Assuming = 0^2 = 

above equation as. 



A/m * —2 


( 5 ) 
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where 


_ ^1 '*'^2 ~*~^3 is the average effective magneton per ion. 


Comparing this equation with asymptotic equation of the curve given by eq.2, one gets. 


—2 


0 = 0^, and Cm = iVppPo? / k 


• 11/2 


or P = ACm / yVpj pO 


( 6 ) 


The experimental value of P can be calculated from the evaluated value of Cm • 
The theoretical values of P,, P^ and Pj are known. Hence, one can obtain the 
theoretical values of P. The experimental and theoretical values of P are given 
in table 3 together with the magnetic ions used to obtain theoretical values of P. 


It is seen from this table that there is a good agreement between the 
theoretical and experimental values of P. This indicates that all the studied 
compounds are essentially ionic and magnetic states of the ions are as indicated 
in the fourth column of the table 3. It can also be noticed that in TbFeCrO^ and 
TbFeMnO^ compounds, there exist CF* and ions and they substitute Fe^* 
ions. This is quite reasonable in view of the natural valency of these elements. 


Table 3. Magnetic ions and average effective magneton number per ion (P) 


Compounds 

Values of P 

Theo. Expt. 

Magnetic ions 

TbFefi, • 

7.14 

7.05 

W^F^*F^^ 

TbFeCrO, 

6.66 

6.61 

r 

TbFeMnO. 

4 

6.88 

6.79 


TbFeCoO, 

4 

6.93 

6.92 






U. P. SINGH, A. N. THAKUR AND V. SINGH 


61 



Fig. 1. Plots of inverse of magnetic susceptibility against absolute ten^erature for TbFej O^, TbFeCrO^, 
TbFeMnO, and TbFeCoO, 
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COINCIDENCES AND FIXED POINTS 
OF MEIR-KEELER TYPE CONTRACTIVE MAPPINGS 
IN RANDOM NORMED SPACES 
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ABSTRACT 

In this paper we introduce the notion of Meir-Keeler type contractive condition for four mappings 
in a random normed space and obtain a coincidence and fixed point theorem for these mappings. Our 
result improves and extends some of the known results in metric, Menger and random normed spaces. 

Key words and Phrases: Random normed space, Menger space, fixed point. 

Mathematics Subject Classifications (2000): 54H25,54E70,47H10,46S50. 

INTRODUCTION 

Serstnev [11] introduced the concept of a random normed space (/?N-space) 
as a probabilistic generalization of the concept of normed linear spaces. The study 
of fixed point theorems in /?A^-spaces was initiated by Bocsan [1] and further results 
have been established in [2,3]. On the other hand, Meir and Keeler [6] presented the 
following remarkable generalization of the celebrated Banach contraction principle: 

Let/ be a self mapping of a metric space (M, d). If for an e > 0, there exists a 
5 > 0 such that e ^ d(u, v) < e + 6 implies d (fu.fv) < e for every u, v in M, then / has 
a fixed point. 

This result has been studied extensively and extended by several authors (see, 
for instance, [5], [7-10]). An excellent bibliography is given by Jachymski [4]. 

The notion of Meir-Keeler type contractive condition for a triplet of mappings 
on a Menger space was introduced in [13] and some coincidence and fixed point 
theorems were proved. In this paper we introduce the notion of Meir-Keeler type 
contractive condition for a quadruplet of mappings on an arbitrary set with values in 
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an RN-space and prove a coincidence theorem for such mappings. We further, derive 
the existence of common fixed point of these mappings commuting only at their 
coincidences. These results, being new in this direction, are interesting generalizations 
and improvement of a number of known results in different settings. 

PRELIMINARIES 

Let £ denote the set of all distribution functions. An RN-space is a triplet 
(5, t) consisting of a real or complex linear space S, the mapping ^ :S L and 
the t- norm t satisfying the following conditions in which ^(u) is denoted by , and 
/^,(.v) is interpreted as the probability of the norm of u being less than jc: 

(a) £;,(0) = 0for all ueS-, 

(b) (;r) = // (x) iff M = 0, where H (a:) = 

(c) = Fuix 1 1?]), where q is a non-zero scalar and x eR; 

(d) f„,,U + }')>/{F„(jc)./=’„(y)) forallM,vin5and;r,y > 0; 

(e) /(x, v) > max {jr + y-1, o} foralljr.y e [0,1]. 

For detailed topological preliminaries of/2A^-spaces, we refer to [11]. 

Suppose the mappings/ g, h and ^ on an arbitrary set with values in an RN- 
space satisfy the following condition : 

Given e > 0, there exists a 8 > 0 such that 

G„.,(e) = 0 and G„_/8 + 6) > 1 - X imply (e) > 1 - , where 

G„..„(e)= min 

^.-*,(2e).F,,.,.(2£)}. (1) 

The condition introduced above will be referred to as Meir-Keeler contractive 
condition for a quadruplet of mappings and it is evidently equivalent to 


J 0 , 

I 1, x>0, 
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Ffy-gyi^) > min{F,„.*,(E),/;i,_*„(E),F^.*,( e), Fj^,^i2e),F^_^i2e )}. (2) 

Throughout this paper X stands for an arbitrary nonempty set and S for an RN- 
space. By C{f. h) we denote the set of coincidence points of/and h, i.e. Of, h) = 
{r ex:fr = hr}. 

The following result for f,g,h,k:X-^S satisfying the most general contractive 
condition of Meir-Keeler type ((1) or (2)) improves, among others, the results of [3] 
and [12] on one hand and generalizes the results of [13] on the other. 

Theorem. Let/ g, A, ^ .• A' -► 5 be a Meir- Keeler type contractive quadruplet such that 
f{X)ugiX)czhiX)nkiX); (3) 

h(X) n KiX) is a complete subspace of S. (4) 

Then 

(i) /and h have a coincidence point; 

(ii) g and h have a coincidence point. 

Further, if X=S,fhu = hju, u ^Cif.h) and gkv = kgv, v€C(g, A) 
then 

(iii) f.g.h and k have a unique common fixed point. 


Proof. For € X, we define the sequences {«„} and {v, }as follows: 

= Vj,. « = 0, 1, 2,... . 

This can be done by virtue of (3). Now from (2), 


= F„ 


(e). 


Similarly 




F , (e) > F„ , (e). 

Thus (s) is a monotone increasing sequence. Since F is a distribution function, 
sup F, (e) = 1. Therefore F, (e) 1. 

To prove that {vj is a Cauchy sequence, suppose the contrary to be true. Then 
there exists s' (= 2e) such that for each integer ^and integers m.n>Nv/e have 
F,,-,.(2e)>l-A., 0<>,<1. (5) 

We choose 8, 0<8 <E,such that (1) is satisfied. Since 
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-v.., (S ^ 6) 1, there exists a positive integer N= N(e,X) 

such that for m, n ^ N 

and/;^.j6+5)<l-X. (6) 

Let (6) be true for an even n. Then (e + S / 3) < 1 - X,, by (5). 

Thus we have a smallest odd integer j>n such that 

F,., (8 + 8/3)<l-XandF,^., je+5/3)>l-X. (7) 

Therefore from (6) and (7), 

F,. (e + 2fi / 3) i min{. (e + 8 / 3).F,,, (8 / 6).F,, (8 / 6)} 

> Tnin{I-X. I-X, 1-X} = 1-X. 

Now G,, , (8 + 6) = , (8 + 8) 

= min {F,(8+8),F, ^^.,^(e+5),F,^ (8+8), 

F,.,.J2(8 + 8)),F,^^..,/2(8 + 8))} 

> min (8+28 / 3),F,^^,.,.^ (8 / 6),F,_^ (8 / 6), 

^v,.-v/8+8),F,^.j8 + 8).F,„^..j8+8),F,,.,/8+8)}>l.;,. 

Therefore from (1), 

F/v (s) = F, (8) = F^ (e) > 1 - X. 

So 

F,,(8 + 8 / 3) ^ min^ (8 / 6),(e),(8 / 6)} > i _ x, 
contradicting (7). 

Therefore {v,} is a cauchy sequence inh(X)n k (A). Let z be the limit. Then there 
exist point r, s such that r = h~'z and $ = k~'z »giving z = hr. 

To prove fr = z, let C/)j.(e,X) be any neighbourhood of fr. Since -> hr, 
there exists an integer N' (e, X) such that 2n ^ N' implies 

> I - and >!->.. (8) 

Now by (2), 
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> (e / / 2). 

Thus hr-z=fr. Similarly ks = z =gs. 

To prove (iii), we note that r e C(f. h), s e Cig, k) and 
(a') hr =fr = z = ks = gs, 

(b') fa =flir = hfr = hz, 

(c’) gz=^gks= kgs = kz. 

Now replacing u and v by r and z respectively in (2) and invoking (a')-(c'), we get 
> ^z-gzi^)y implying, gz = z. By(c'), 

kz = gz = z. 

Similarly taking u = z, v = s,we obtain fa=z. 

Thus by (b') and (o') 

fa - hz - z = gz = kz. 

The uniqueness of z as the common fixed point off, g.h and k can be easily seen 
from (2). 

Taking X=S,h=k = /and S to be orbitally complete (see,[12]) in the above 
result, we get the following corollary. 

Corollary. Let the self mappings f, g on an (f, g)- orbitally complete RN-space S 
satisfy 

^ju-gyX^) > min{F„.,(e),F^.„(8),F^.,(8), F,„_,(28),F^,.„(28)}. 

Then for every i^eS, either 

/or g has a fixed point in the (f, g)- orbit of or 

/and g have a unique common fixed point z such that lim u„ = z. 

The above result for a pair of m^pings is an extension to F^-spaces of Theorem 
2 of [10] under the most general contractive condition of Meir-Keeler type. This 
result is also an improvement over some of the known results, including that of [10], 
in the sense that we do not use the continuity of either mapping here. 
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ABSTRACT 

In this paper we prove some fixed point theorems for densifying mappings on probabilistic 
metric spaces. Several results for densifying mappings in metric and probabilistic metric spaces 
are obtained as special cases. 

Mathematics Subjects Classifications (2000): 54H2S, S4E70,47H10 

Key-words and Phrases: Probabilistic densifying mappings, probabilistic weakly 
(|) -contractive, fixed point. 

INTRODUCTION 

The notion of densifying (also called condensing) mappings was introduced by 
Sadovskii [11]. Furi and Vingoli [4] studied the fundamental properties of Kuratowski 
measure of noncompactness of a bounded set in a metric space. Later, several 
mathematicians investigatedfixed point theorems for densifying mappings [6,8,9, 
10,13, 14]. 

Bocsan and Constantin [3] extended Kuratowski measure of non-compactness 
of probabilistic metric spaces (PM-spaces). The concept of probabilistic densifying 
mapping was introduced by Bocsan [ 1 ] (cf [2]), and some fixed poing theorems for 
these mappings have been proved [2,5,7,12,13]. 

In this paper we introduce the notion of (p; q^, q^, ., qj- probabilistic 

densifying mapping and prove a fixed point theorem for such mappings. Further we 
extend this result for a pair of mappings. The paper is concluded with a result which is 
the generalization of a result of Pant, Tivari and Singh [7]. Several fixed point theorems 
for densifying mappings on metric and PM-spaces may be obtained as special cases 
of our results. 
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DEFINITIONS AND RESULTS 

Definition 1 [ 1 ]. Let ^4 be a nonempty subset of A', A function (.) defined by 

is called the probabilistic diameter of A. 

A is said to be bounded if sup D^{x) = 1, R being the set of real numbers. 

x<R 

Definition 2 [3]. For a probabilistic bounded subset .,4 ofA', the function a^(.) defined 

bya^(jc) = supje >013 a finite cover A of A such that D,(X)^e for all 5 e A} is 
called Kuratowski function. 

The following properties of are proved in [10]. 

(i) for every probabilistic bounded subset A of A'and for all ^ e /? J 

(ii) OL^(X) > D^(X) for every probabilistic bounded subset A of A' and for all 
X eRi 

(iii) If (|) ^ A c c A then a^(x) > ag(x) for all x eR, 

(iv) If A and B are probabilistic bounded, then for each x eR 

^ rain {a^(A),a(A)}; 

(v) If A is a probabilistic bounded subset of X, then 

a 

A 

(vi) A probabilistic bounded subset A of X is probabilistic precompact iff 
a_ =//; 

A 

(vii) If the t-norm t is 't-min' then ® ® COA for every probabilistic bounded 

subset A of X, where denotes the closed convex hull of A in the (e,X) 

-topology on X. 

Definition 3 [ 1 ]. Let (X, F) be a probabilistic metric space. A continuous mapping f of 
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X into X is called a probabilistic densifying mapping iff for every subset A ofX such 
that 

< // we have >a^ 

We introduce the following: 

Definition 4. A mapping/: X-> is said to be .«»,) 

probabilistic densilying if for >1 c X, 

(1) /'’ is continuous and 

(2) I o, 

/«l 

whenever^flyOt^, .^^^ <H where p,q^ q^...q^ are all non-negative 

m 

integers and the a- s are non- negative reals such that X ~ 

/«1 

Definition 5. A mapping/: A"is called probabilistic weakly 

0 _ contractive iff 

(3) for every/4c.If,<//implies a < a ^ 

(4) for every u, v eX. u:^v 
0 (/(w),y(v))>0 (m.v) 

where 0 is a x-Continuous mapping [3] of XxX into It is to be noted that 
T - continuity coincides with upper semicontinuity [5]. 

THEOREM 1. Let /be a (p; ^j... - probabilistic densifying mapping defined 

on a complete Menger space (X, I, t) wherein t (x, a:) = 1 such that /'’is 

probabilistic weakly 0-contractive. If, for some eX, the 

sequence of iterates {m„ }is bounded, then/has a unique fixed point inX. 

PROOF. Let/4 = U{m„} where u„ = fu„,,n = 1,2,....then 

»=o " 

r{A)=f'\_r\ . /'''{^) = /’■«.../’"■««. 

Thus f{A) and fj {A),j = 1, 2,..., m all differ from A only by a finite number of 
terms. Hence 
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If I thenby(2), 

m 

a/’'’*'*'> y a a =a „ which is, a contradiction. 

i—J J (A) i' (A) 

h\ 

Thereforeis compact. 

We consider vj/ : f'‘{A)^C defined by 

M/(h) = 0(«,/'■(«)). 

Then vj/ is u.s.c and attains a maximum value at some point, s2iy,zef'’{A). 

The continuity of f gives 

C /'■(/'’(A)) = fP(A) c i.e. r(z)G/"(A). 

IfZ?!: /' (z), then 

vj/(/''(z)) = 0(//(z),/''’(z))> 0(z,/'’(z))=/(z). This contradicts the 
maximality of v|; atz. Hence z = /''(z). 

Using the probabilistic weak 0 -contractivity of f" it is easy to see that z is the unique 
fixed point of f’’. 

Also f{x)=f {/’’{z)). Hence z = f{z) by the uniqueness of z. 

REMARK 1. If f is a (1; 0) probabilistic densifying moping, then the above Theorem 
will reduce to the following; 

COROLLARY. Let /be probabilistic densifying and probabilistic weakly 

0-contractive mapping defined on a complete Menger space with 

t{x, jc) = 1. If for some € X, the sequence of iterates {u„ } is bounded, then/ has 

a unique fixed point in X. 

Now we extend Definition 4 for a pair of mappings. 

Definition 6. A pair of mappings 5= {/,g} defined on a Menger space (X. t) with 

t{x, x)= 1 is called . ,q„)-probabilistic 

densifying if for A c X, 

(5) /'■ and g" are continuous and 
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( 6 ) 




> 


y-i 


whenever <W. where 

the a/f are non-negative reals such that 


.are are all non-negative integers and 

m 

= j- 

i-i 


THEOREM 2. Let 5 = {f, g} be a pair of commutative and (p;q,,qi . q„)- 

probabilistic densifying mapping on a complete Menger space (X. t) with 

t (;c, ;c) = 1 such that f<‘g’> be probabilistic weakly 0-contractive. Further, let for 

some ifnisoddand, ifn is even, then/and g have unique 

common fixed point. 

oo 

PROOF. Let=UK}- 

n=0 

Now fi'g'i (A) and (fg)''-(A),j = !,2 .n, all differ from A only finite number of 

terms. Hence, 




m 

If Z (/*)’> (A)^ ’ then by (6) 

y*i 


^ Z “ aZ “ “a ,whichisabsurd. 

Therefore ^ is compact. 

Define v(/: A -» /? such that V|/ (m) = 0 (m, (/Y')(«))• Then vi; is u.s.c. on 
as 0 is u.s.c. So V)/ has a maximal value at z e A . 

Now /'’g*(A) c /'’(g*(A))c A. • Hence /'’^’(z) € A 






74 


FIXED POINTS OF PROBABILISTIC DENSIFYIN6 MAPPINGS 


Ifz ^ f’g-iz) € A, then v (/'VW) = 0(/V/U)./V/V(z)) 

> 0{z,f''g'*(z)) = M/(z), contradiction to themaximality of/atz. Hence /'’g»(z))=z. 

Now we prove the uniqueness of z. Let w be another fixed point of fg ", then 
0(z,h ) = 0 (/''g* 2 ,/'>’w)> 0 ( 2 ,^) which is not possible. 

Therefore z = w. 

Finally, we show that z is the unique common fixed point of f and g as well. 

Notice that z = fg^z) implies f {z)= f(j'’g'>(z)) = f'g^fiz)). 

Thus f (z) is the fixed point of . But z has been shown the unique fixed point of 

/’’g’’ and so z =/(z). 

Similarly z=g (z). Hence z is the unique common fixed point of f and g. 

Definition 7. Let 5 ={/'’, g’}, for some non-negative integer p and q, be a pair of 
self-mappings on Menger space {X.^, t) with / (x, x) =1. Fore ^ the sequence 

L («o )= {Wo./'’“o > } is called the joint sequence of iterates of S' 

at Mo. 

The following lemma is obvious (see also [7, Lemma, p. 152]). 

LEMMA. Let (X. ^,t)\)edi complete Menger space withS^'^P^('*»-*) “ i and / g be 

self mappings of X. Then/and g are probabilistic densifying mappings if /"for every 
pair of probabilistic bounded subset ..4 and B of A', 

whenever 

The following result is an extension of Pant, Tivari and Singh [7, Theorem, p. 

152]. 

THEOREM 3. Let S={} be a pair of commutative probabilistic densifying 
mappings defined on a complete Menger space {X, t) such that supt(x,x) = 1^ 
where p and q are some positive integers. Furhter, let fg^ be probabilistic weakly 0 
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contractive and for some € X, the joint sequence of iterates of S at Mq is 

bounded, then/and g have a unique common fixed point in X. 

PROOF. Let A = {., 

B = .} and C = .}, 

so that A = B u C and A = u {mq}- 

Now replacing f and g by f and g" respectively in the proof of [7, Theorem, p. 

152], we get f'‘g‘> (z) = z. Uniqeness of z can be seen easily. 

Rest of the proof follows from Theorem 2. 

REMARK 2. By taking = 1, ^ = 1 in the above theorem we get the result in [7, 
Theorem, p. 152] as special case. 
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ABSTRACT 

In this paper, we introduce a common fixed point principle which includes several known 
contractive definitions as particular cases and employs a Lipschitz type analogue of known 
contractive definitions. 
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INTRODUCTION 

The study of common fixed point of mappings satisfying contractive type 
conditions has attracted a great deal of research activity during the last two decades. 
The most general of the common fixed point theorems pertain to four mappings, say 
A, B, S and T of a metric space (X, d), and use either a Banach type contractive 
condition of the form 

d{Ax,By)<hm{x,y),0<h<\, (1) 

where m{x,y) = mdx\d{SxJy),d[Ax,Sx),d{By,Ty),\d[Sx,By) + d[Ax,Ty'f^l2^, 

or, a Meir-Keeler type (e,5 )-contractive condition of the form 

given 8 > 0 there exists a 6 > O such that e ^ m[xyy) < e + 6 d{Ax,By) < e, (2) 

or, a (|> -contractive condition of the form d{Ax, By) ^ <|>(w(x,y)), (3) 

involving a contractive gauge function is such that for each t > 0. 

Clearly, condition (1) is a special case of both conditions (2) and (3). A 
(() -contractive condition (3) does not guarantee the existence of a fixed point unless 
some additional condition is assumed. Therefore, to ensure the existence of common 
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fixed point under the contractive condition (3), the following conditions on the function 
4 > have been introduced and used by various authors. 

(I) (|) (/) is non decreasing and //(r- ^ (0) is non increasing [2], 

(II) (|i (0 is non decreasing and lim^ (|> "/(0=0 for each / > 0 [4, 8 ], 

(III) <|) is upper semi continuous [ 1,4,7,11, 1 2] or equivalently, 

(IV) <|) is non decreasing and continuous from right ([ 17]). 

It is now known (e.g, [4,15]) that if any of the conditions (I), (II), (III) or (IV) 
is assumed on <|), then a <j) -contactive conditions (3) implies an analogous (e, 8 )- 
contractive condition ( 2 ) and both the contractive conditions hold simultaneously. 
Similarly, a Meir-Keeler type ( 8 , 8 )-contractive condition does not ensure the 

existence of a fixed point. The following example illustrates that an (e, 8 ) -contractive 
condition of type ( 2 ) neither ensures the existence of a fixed point nor implies an 
analogous (j) -contractive condition (3). 

EXAMPLE 1. ([15]) Let A"= [0,2] and dhe the Euclidean metric coX. Define 

f\X-^Xhyfx = {\+x)l2 \ix<\ \ fx = a ifx>\. Then, it satisfies contractive condition 

8 >max {d{x.y), d{x,fx), d {y,fy), [d{xjy) + d{yjx)]l2] <e + 6 => {fx,fy) < e, 

with 5 (e) = 1 for E > 1 and 5 (e) =1 - e for e > 1 but/does not have a fixed point. 
Also does not satisfy the contractive condition. 

d ifx. fy)<^ (max {d {x, y), d (x, fx), d (y, fy), [d (x, fy)+d (y, fx)] 2 }). 
since the desired function ^ (/) cannot be defined at / = 1 . 

Hence, the two type of contractive conditions (2) and (3) are independent of 
each other. Thus, to ensure the existence of common fixed point under the contractive 
condition ( 2 ), the following condition on the function $ have been introduced and 
used by various authors. 

(V) 5 is non decreasing [ 10,11 ] 

(VI) 5 is lower semi continuous [5, 6 ]. 
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Jachymski [4] has shown that the (e,5)-contractive condition (2) with a non 
decreasing § implies a (|> -contractive condition (3). Also, Pant etal.[\S\ have shown 
that the (£,8)-contractive condition (2) with a lower semi continuous 6, implies a 
(|) -contractive condition (3). Thus, we see that if additional conditions are assumed 
on 8 then the (e,8) contractive condition (2) implies an analogous (|> -contractive 
condition (3) and both the contractive conditions hold simultaneously. 

It is thus clear that contractive conditions (2) and (3) hold simultaneously 
whenever (2) or (3) is assumed with additional condition on 8 or(|> respectively. It 
follows, therefore, that the known common fixed point theorems can be extended and 
generalized if instead of assuming one of the contractive condition (2) or (3) with 

additional conditions on 8 and ^, we assume contractive condition (2) together with 
a Lipschitz type analogue of condition (3); that is, a condition of the form 

d{Ax, By) < max|fc,^i(5jc, Ty),k 2 ^d{Ax, Sx) + d[By, Ty)] / 2,^d{Sx,By) + d{Axy Ty)]/ 2| 

it, > 0 and 1 < A '2 < 2 . We prove a common fixed point theorem for four mappings 
using this approach. It is independent of known earlier results. It may be noted that all 
the known results have been dealt with the case /:, = !; ^ ^o all such 

results are obtained as special Case of our theorem. 

Two self-mappings A and 5 of a metric space {X, d) are called compatible (see 
Jungck[5])if Lim d(ASx„,SAx„) = 0, whenever {xj is a sequence in^such that 

/I 

Lim Ax„ = LimSx„ =t for some tin It is easy to see that compatible maps commute 
at their coincidence points. 

RESULTS 

We prove the following theorem with the notation M(x,y) defined as 
M[x,y) = ma\^d{Sx,Ty),d{Ax,Sx),d{By,Ty),[d{Sx,By) + d{Ax,Ty)y2) 

Theorem 1. Let (A, S) and (B, T) be compatible pairs of self mappings of a complete 
metric space (X, d) such that 


(i) AXclTX,BXczSX, 
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(ii) given e > 0 tiiereexists 6 > 0 such that e < M{x,y) <e + 6=> d{Ax,By) < e, and 

(iii) d[Ax,By) < mdix.^k^d{Sx,Ty),k^^d[Ax,Sx) + d{By,Ty)^l2 

\d{Sx,By) + d{Ax,Ty)\l2}, for >0,1<^2 <2. 

If one of the mappings A, B. S and T is continuous then A, B, S and Thave 
unique common fixed point. 

Proof. Let x^ be any pointing Define sequences and {y„} inA'givenbythe 
nile y 2 „ = Ax 2 „ = (4) 

This can be done by virtue of (i). Then following the argument by Rhoades et al. [ 19] 
it follows that {yj is a Cauchy sequence inX But A^is complete so there exists a 
point z in A' such that y„ -> z. Also, using (4), we have 

yzn = ^ “ ^21*2 (5) 

Suppose that S is continuous. Then SSx 2 „ -> Sz, SAx 2 „ -> Sz and compatibility 
of A and S implies that ASx 2 „ Sz . Also, since AXc: TX, corresponding to each 
value ofn, there exists, Zj, inA'suchthat/IS'Xj, = Tz^,. Thus^45x2,= Tzj^-^iSz and 
SSx 2 „ Sz. We show that Hm, = Sz. If not, then there exists a 
subsequence {Bzj„} of [Bz 2 „) a number r >0 and a positive integer Wsuch that for 
each m ^ N ,weha.ved (ASx 2 „ Bz 2 „) ^ n d{Sz,Bz 2 „)'t randinviewof(iii), weget 

</(A5jc 2 fiZ,,„) > max [k^d{SSx 2 „JZi„\k 2 [diASx 2 „,SSx 2 „ )+ diBz 2 „,Tz 2„)1 /2, 

[diSSx2„_Bz2 „) + diASx2 „, Tz2 „)] / 2 }, 

which, on letting m oo, yields 

diSz, Bz 2 „ ) < [d(Sz, Bz 2 „ )] / 2 < d(Sz, Bz 2 „), a contradiction. 

Hence lim„_,, Bz 2 „ = iz. We claim that i4z=5'z. If i4z9i&, then by virtue of (iii), for 
sufficiently large values of n, we get 

d(Az,Bz2„) < max[k^d(Sz,Tz2„),k2[diAz,Sz) +d{Bz2jz2n)]/'^^ 

[d(.Sz,Bz2j + d(Az,Tz2j]/2}. 




R. P. PANT AND K. JHA 


81 


On letting ► oo, this yields d (Az, Sz) < \d (Az. Sz)]/2 < d (Az, Sz), a contradiction. 
Hence Az = Sz. Also, since AX c TX, there exists a point w in such that Az = Tvv. 
If 7W, using (iii) we get 

d{Az, Bw) < max {k^d{Sz, Tw), [d{A 2 , Sz) + d(Bw, TV)] / 2, 

[d{Sz, Bw)+d{Az. 7V)]/2}, 

= [d (Bw, TV)]/ 2<d (Bw, Tw) = d (Bw, Az), a contradiction. 

Hence Az = Bw and so, Sz=Az = Tw- Bw. 

Since compatible maps commute at their coincidence points, we get ASz = SAz 
and BTw = TBw. Moreover, AAz = ASz - SAz = SSz and BBw = BTw = TBw = 
TTw. If Az i^AAz, using (ii), we find 

d (Az, AAz) = d (AAz, Bw) < max {d (SAz, Tw), d (AAz, SAz), d (Bw, Tw), 

[d (SAz, Bw) + d (AAz, TV)]/2}, 

= d (Bw, AAz) < d (AAz, Bw), a contradiction. 

So that Az = AAz = SAz and so Az is a common fixed point of ^ and S. Similarly, 

Bw (= Az) is a common fixed point of B and T. Uniqueness of the common fixed point 
follows from (ii). The proof is similar when T is assumed continuous in place of S. 
Moreover, since AX c TX and BX c SX, the proof follows on similar lines when A 
or 5 is assumed to be continuous. This establishes the theorem. 

We now give an example to illustrate the above theorem. 

EXAMPLE2. LetX= [2,20] and d be the Euclidean metric onX Define yl, B, S and 

T. X-^XdiS follows: 

Ax = 2 for each x 

5x = 2ifx<4or>5, Bx = 3+xif4<x<5; 

Sx = xifx<8, Sx = 8 ifx>8; 

Tx = 2, ifx< 4 or > 5, Tx = 9+ x if4 ^ x< 5. 

Then A, B, S and T satisfy all the conditions of the above theorem and have a 
unique common fixed point x = 2. It can be seen in this example that A, B, S and T 
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satisfy the condition (ii) when 6 (e) = 1 if e > 6 and 5(8) = 6 - eif 8 ^ 6 . Thus, 6 (e) 
is neither non decreasing nor lower semi continuous. It can also be verified that the 
mappings A, B, S and T do satisfy the contractive condition (iii) with k^ = l and 
/tj = 1 . However, A, B. S, and T do not satisfy the <{> -contractive condition (3) since 
the required function <|> (/) can not be defined at / = 6 . Hence we see that the present 
example does not satisfy the conditions of any previously known common fixed point 
theorem for contractive type mappings, since neither the mappings satisfy a 
(|> -contractive condition nor 6 is lower semi continuous or non decreasing. 

REMARKS. As various assumptions either on <|> or on 6 have been considered to 
ensure the existence of common fixed points under contractive conditions, so our 
Theorem 1 improves the results of Boyd and Wong [1]. Carbone et. al. [2], 
Matkowski, [ 8 ], Pant [10, 11, 12], Pant and Pant [14], Park and Rhoades [17], 
Singh and Kashahara [20], Jungck [5], Jungck et al. [ 6 ]. Jachymski [4], Maiti and 
Pal [7], and Park and Bae [16] for the case when k^=\ and k = 1 . All such results are 
obtained as special case of Theorem 1 when ^,=1 and = 1. Also our theorem, 
thus, generalizes all other similar results for fixed points and allows k, to take values 
other than 1 by taking a Lipschitz type contractive condition. 
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ABSTRACT 

The object of this paper is to develop a non-equilibrium thermodynamic model of a 
population ecosystem and to study the thermodynamic criteria of stability of the ecosystem. 
Mathematics Subject Classification (2000) : 92Ai7 

Key words: Ecosystems, Non-equilibrium thermodynamic model, Lyapunov function, entropy 
function, ecological stability. 

INTRODUCTION 

An ecosystem consisting of an enormous number of interacting elements in 
relation to its environment is a complex system [1,2,3]. The statistical mechanical 
model of ecosystems initiated by Kemer [4,5] based on Gibbsian ensemble theory, 
although very powerful, has however some limitations. It is limited to conservative 
system, but most of the ecological model equations are of non-conservative type. In 
recent times great efforts are going on in developing general non-equilibrium 
thermodynamic models of complex ecosystems [6,7,8,9,10,11,12,13,14,15, 
16]. The thermodynamics has the ability to provide a unified description of macroscopic 
systems which are independent of the details of ongoing processes in terms of state 
functionals depending on limited number of observables. 

In an earlier paper [7] we have developed a non-equilibrium thennodynamic model 
of Lokta-Volterra ecosystems. The object of the present paper is to supplement the earlier 
one. In this objective our first task is to derive the expression for entropy-production of a 
general model ecosystem, which is independent of Lokta-Volterra model equations or any 
other model equations ofbiologjcal and chemical reaction systems. The second is to explore 
the extent to which the thermodynamic model can provide useful characterization of the 
dynamical complexity of stability of some model ecosystems. 

ECOSYSTEM: THERMODYNAMIC MODEL 

In applying thermodynamics to the theory of ecosystem we assume that biological 
phenomena, although not completely determined by the laws of macroscopic physics 
and chemistry, are not in contradiction to them [14]. For a consistent thermodynamic 
theory of ecosystems, we make the following assumptions: 

1. We assume the ecosystem as a regional community of biological species 
(excluding biotic components) embedded in an environment. 

♦ Department of Applied Mathematics, University of Calcutta, Kolkata-700 009. India 
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2. Let /V,, (/■ = 1 ,2,.be the population of the i-th species and N, is assumed 
to he governed by the system of Kolmogorov's type of equations: 

dN 

.(/=1.2. n) (2.1) 

where the functions Fi's are assumed to be continuous function of 

having partial derivatives in some open set Q = [iV, ^ 0 : j = 1,2.«]. 

3. We assume that the energy of the ecosystem is stored exclusively within the 
organisms so that the transfer or flow of energy is governed exclusively by the system 
of dynamical equations characterizing the model ecosystem. 

4. The model equations are taken as the basic equations for the non-equilibrium 
thermodynamic modelling of the ecosystem. This is analogous to Kemer's statistical 
mechanical modelling of ecosystem on the basis of Lotka-Volterra equations. 

5. The state of ecosystem is assumed to be defined by the set of population sizes 

N = [N,,yV 2 ,...,Njand they are assumed to be extensive variables for the 
thermodynamic description of the system. 

Let N* = be the stationary values of N 

Let S be the entropy of the system at any time t and is assumed to be a function of the 
state (or extensive) variables , that is, S = S{N^,N 2 ,..-N„). Let 

AiS be the change of entropy of the system in transition from the stationary state 
N* = , A ^2 ,..., A^* j to the non-stationary state N = [//,, A ^2 ,..., ], that is, 

AS = S-S„ (2.2) 

Let us now develop a non-linear thermodynamic model of the ecosystem 
described by the dynamical equation (2.1). Expanding about 

the stationary state [ A^,/V^,....iV*] we can reduce the system of equation (2.1) 
to the form 

- (Nj-N]) 

jL 


1 dN, 
M dt 


=/,(A^;,A^;.A^;)+Z 




M 

dN 


1 /I n 

4 x 5 : 

^ 7*1 km\ 




dNjdN, 


(Nj-N]){N,-N:y 


(2.3) 
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Neglecting higher-order terms of bNj = Nj - N], we have 


1 dN, ^ 
dt ^ 


M. 

dNj 


K-^;) 


Let us now define the thermodynamic fluxes the specific growth rate 


(2.4) 


J, = 


1 dN. 


N, dt dt 


(In A/,).(i= 1.2.n) 


(2.5) 


and the corresponding thermodynamic forces Xj as 

X^ = N^-Nl{i = \,2,..,n) (2.6) 

which are the deviation of the macroscopic or state-variables N. from their stationary 
value Ni . With these definitions of thermodynamic fluxes, y, and forces Xj the system 
of equations (2.4) reduces to the forni of linear phenomenological relations of irreversible 
thermodynamics [10]: 


. «) 

M 


(2.7) 


where 




Ml 

aw, 




( 2 . 8 ) 


are the phenomenological coefficients. We now find the explict form of the 
entropy-production AS . According to Onsagar the entropy-production is given by 
[ 10 ] 


dt 





/«! 


(2.9) 


Integrating (2.9) from the stationary state [W,*, Wj.. Wjj,] to the non-stationaiy state 

[W, Wj.WJ, we have 


AS = ± 


(N,-N;)-<-N;inL 


( 2 . 10 ) 


which is the expression of the change of entropy (except a multiplicative constant) in 
transition from the stationary state [Wj‘,W^,....,W|[] to the non-stationary state 
[W,, ...., WJ.Theexptession(2.10)playsaCTUcialiDleinlhestiKfyofstabililyofthesystem. 
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ENTROPY AND LYAPUNOV FUNCTION; 
THERMODYNAMIC CRITERIA OF STABILITY 


A class of Lyapunow function for the dynamical system (2,1) is given by [ 17] 


i=SA/;’0 


''w'' 


/*! 




(3.1) 


where 0 is a convex function of N, having isolated minimum at Nj-Nl and increasing 
monotonically. Now by Lyapunov theorem [17] the stationary state [] 
is stable if dL/dt > 0 taken along the trajectory of (2.1) and unstable if dU dt> 0 
even for one trajectory of (2.1). If we consider the function 0 in (3.1) by 


P -IL 

) - In - 1 , where s, - 

Then the corresponding Lyapunov function will be given by 


(3.2) 




{N,-N])In 


K. 

N, 


(3.3) 


which is the same expression of entropy change (2.10) derived fix)m the non-equilibrium 
thermodynamic modelling of the dynamical system (2.1). The entropy change 

given by (2.10) which is positive and vanishes identically for Ni = N*, (i = 1,2,...,«) is 

this a Lyapunov function defined over the state space [ A^,, A/j,..., ] of the system. 

For the system of model equations (2.1) and the entropy change a5 as the Lyapunov 
we have the following theorem: 

Theorem: The model system desribed by the system of equations (2.1) is globally 
stable in the positive orthant of the state space [ V,, ,..., N, ] if the time derivative 
along the trajectory 


^(A5) = (3,4) 

i=I 

in the positive orthant and does not vanish identically along any sholution of the system 
( 2 . 1 ) other than the stationary state [iVf, ,..., ^* ]. 


The proof follows directly from LaSalle's extension of the direct. Lyapunov 
theorem. The equality in (3.4) corresponds to the marginal stable state having 
indeterminate stability characteristing stable from the unstable regions [17]. 
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Let us now derive the thermodynamic criteria of stability from the above general 
criteria (3.4). Applying Taylor’s series expansion of /(A/,, ,,.., Aj and using (2.8) 
and (2.9), we have 


£ 

dt 


(a5)=2ZK-^;)| 

/=! k=\ 




, /•=i *=i 


(3.5) 


The criteria of stability of the stationary state [at* , AT*,..., A’ ] then becomes 

the negative-definiteness of the quadratic form (3.5) or of the matrix of the 

phenomenological coefficients. In the next section we shall test this criterion of stability 
for some model ecosystems. 

Let usnow consider the stability of Lotka-Volterra ecosystem in the light of the 
above criteria of stability. We consider an ecosystem governed by the generalized 
Lotka-Volterra system of equation [7] 


dN, 

dt 




(i = 1,2.n) 


(3.6) 


where a,, (/ = 1,2,..., n) may have any sign and there is no restriction on the coefficients 

b,j,{i,j = 1,2,..., n). The non-trivial stationary state [aT* , Vj,..., iV’ ] of the system 
(3.6) is given by 


'Zb„N]=a,, (/ = 1,2 . n) 

/ = ! 


(3.7) 


The necessary and sufficient condition for the uniqu solution of (3.7) with all Nj 

positive is that the number of species n must be even and that all a\s be not of the 
same sign [18]. With the help of (3.7) the system of equation (3.6) can be written as 


1 dNt 

N, dt 


/*=! 


(3.8) 


The right hand side of (3.8) is already a linear function of [VpV,*-••»”]• 
Comparing with (2.7), we see that 

(/,./ = l,2,...,n) (3.9) 

so that the criteria of stability of the stationary state [n* for the generalized 
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Lotka-Voltera system (3.6) is the positive definiteness of the matrix [6^] [18]. The 
Lotka-Volterra system (3.6) written in the form (2.1) is a special form of the system, 
the functions in the right hand side of(2.1) are simply linear function 

of the state variable. But all ecological systems are not of Lotka-Volterra type and the 
functions /(Vp A/j,..., A/,) are in general non-linear. In the next section we shall 

consider these two types of systems namely Lotka-Volterra type and non-Lotka- 
Volterra type to test the thermodynamic criteria of stability. 


APPLICATIONS: STABILITY OF SOME MODEL ECOSYSTEMS 

In this section we shall study the problems of stability of some model ecosystems 
on the basis of the theories developed in the previous section. 

A. Competitive system of two species: 

let us first consider a two species Lotka-Voltka competitive model ecosystem 
described by the system of equations 

i^=N,F,{N„N,) (4.1) 

We linearize the system of equations (4.1) about he stationary state (V,* ,Nl), 
Proceeding as before, we have the system of phenomenological equations 

y, = J, = 4.^, +^2^2 (4-2) 

where 


The criteria of stability of the stationary state (V,*, N ]) is given by, L, i < 0, £,,42 “ 
4:41 > 0. This last condition can be writtQi in terms of partial derivatives of ^(Wp.V 2 ) 
and 4 (^p ^ 2 ) respect to V, and evaluated at [N\ , Vj ) as follows; 




[ 54 ] 

> 

[54] 


[54] 

_av,_ 

St 


St 

dN2_ 

St 



(4.4) 


This is the criterion of conq)etetive exclusion introduced first by Gilpin and Justice 
[ 19] without any proof Let us now consider in particular a competetive ecosystem of 
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the two species governed by the system of equations 

^=N,{a-bN,-cN,) (4,5a) 

dN 

-^=N,{e-JN,-gN,) (4,5b) 


where a,b,c,e,f,g are positive constants. The stationary states (iV*, ) are given by, 

(i) = (0,0), (both the species die out). 

(ii) [N\,NI)~ (0,e/g), (first species die out). 

(iii) ) = (a/b,0), (second species die out). 


ag~ec af -be 


ygb-cf' fc-bg 


, (both species exit). 


defining the thermodynamic fluxes J. and forces Xf as in (2.5) and (2.6), the 
system of equations (4.5) can be reduced to the form of linear phenomenological 
relations; 

y, = *^2 ~ Al-^l ^2-^2 

where 


J - 

\dF,(N„N,)' 

^11 - 

dN, 

r — 


^12 - 

dN, 

4.= 


dN, 

L » Jj 

42 = 

'dF,{N,.N,)'] 

dN, 


(4.6a) 


(4.6b) 


(4.6c) 


(4.6d) 
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In this case are independent of different values of Here 

Z,,, = < 0 as 'b' is positive. Now, Z.,,Z .22 ~ > 0, that is of bg > cf. 


B. A Prey-Predator System: 

As a second illustative example let us consider a non-Lotka-Volterra of model 
prey-predator system described by the system of equations [20,21 ] 

jf^^=Nfi-N,)-N,=F,{N„N,) (4.7a) 

= (4.7b) 


where A^, and N 2 are the populations ofthe prey and predator species respectively. 
The stationary states (A^*, A^ 2 ) given by 

(i) (A^i* ,Nl)- (0,0), (both prey and predator die out). 

(ii) (z/*,A^*) =(1,0),(onlypredatorspeciesdieout). 

(iii) (A^*,A^*)= (p,(i-p^), (both prey and predator species coexist). 

Defining the thermodynamic fluxes y, and fores X. as in (2.5) and (2.6), the 
system of equations (4.7) can be reduced to the form of linear phenomenological 
relations 


J^ — + L^2■X2f J2 ~ ^22^2 


where 



=\-2n: 


= -c 


(4.8a) 


(4.8b) 


jt/ 


4.= 


' 5F2(n„N2) 

dN, 


= -/ 


(4.8c) 



C.G. CHAKRABARTI, MALAY BANOYOPADHYAY AND RAKHi BHAHACHARYA 


93 


42 = 


dF,(N,,N,) 

dN. 


= -g 


(4.8d) 


Let us now test the criteria of stability for different steady- states. The steady- 
state (0,0) is trivial and unstable. We consider the second steady-state (1,0), for this 

the phenomenological coefficients 4/ are given by = -I,/-,: = ~J-4i = = 0. 

Here L,, <0,4ii42 “ 424i • So the steady-state (1,0) is stable ^> 0 -Letus 


now test the steady-state (p, p - ). For this interior steady-state L.. is given by 

41 “ ^ ~ 42 ~ 4l = ^22 ~ ^ 

So the steady-state (p,p - p‘) is stable if L,, = 1 - 2p < 0 or p < 0.5 and 4i42 “ 


Z,, ,Z- 2 , =k>0. Thus for ^ > 0 the steady-state (p, p - p^) is a stable for p > Oi and 

unstable if p > 0.5. For p = 0.5, ^ = I the state is marginally stable. In fact, the parametric 
value p = 0.5 is a Hopf-bifurcation point [21,22]. 


CONCLUSION 

The concept of stability is very fiindamental in the study of biological communities 
and there are a lot of theories on stability for ecological communities [18,23,24]. A 
community is considered stable when the component populations do not undergo 
sharp fluctuation and this definition is closer to the thermodynamic (or more correctly 
to the statistical physics) notion of system stability. Evidently general thermodynamic 
concepts (for instance the second law of thermodynamics for closed syste and its 
extension for open systems) should be applicable to biological and in particular, to 
ecological system. In fact, the different thermodynamic quantities like entropy, 
thermodynamic potential etc. are Lyapubnov function and the stability analysis based 
on these thermodynamic quantities are special cases of general dynamical analysis of 
stability based on Lyapunov function [25]. The present paper is an attempt in this 
direction. On the basis of the thermodynamic model of a general ecological system 
we have derived the thermodynamic criteria of ecological stability and tested these 
criteria for some model ecosystems. 
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PERIODATE OXIDATION OF O-CHLOROANILINE - A 
KINETIC MECHANISTIC STUDY 

R. D. KAUSHIK*, A. K. CHAUBEY* AND R. K. ARYA* 

(Received 20.12.2002) 

ABSTRACT 

The kinetic-mechanistic studies related to the periodate oxidation of 
o-chloroaniline in acetone-water medium have been reported. The studies were made by 
employing the spectrophotometric technique. The order w.r.t. both oxidant and substrate 
has been found to be one in each. The rate increases on increasing the ionic strength and 
decreases with decrease in dielectric constant. The rate-pH profile has been discussed. 
The effect of free radical scavengers rejects the possibility of a free radical mechanism. 
The thermodynamic parameters have been evaluated and discussed. The main product of 
oxidation characterized by melting point and UV-VIS, I.R. and N.M.R. spectrum was 
chloro-1,4 benzoquinone. A suitable mechanism has been proposed and the rate law 
derived. 1 mol of the substrate was found to react with 2 moles of the oxidant. 

Key word & phrases: Periodate oxidation ; o-chloroaniline; kinetics and mechanism. 

INTRODUCTION 

The kinetic studies made on the non-Malapradian oxidation of aromatic 
amines by periodate ion are rather few [4,7,11,12,14]. In continuation to 
earlier studies [8,9], the resutls of the periodate oxidation of o-chloroaniline 
(OCA) are being reported and discussed. 

MATERIALS AND METHODS 

OCA and sodium metaperiodate of E.Merck A.R. grade were used after 
Zn-dust distillation/recrystallization. Doubly distilled water and other 
chemicals of CDH A.R. grade were used. Thiel, Schultz and Koch buffer 
[2] was used for maintaining the pH of the reaction mixtures. 
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The progress of the reaction was followed by recording the absorbances 
of the light orange colored reaction mixture on Schimadzu double beam 
spectrophotometer, UV-150-02 at the wavelength 460 nm i.e. the of 
the reaction mixture, keeping the pH at 5.0 (except for the kinetic runs in 
which the effect of pH was being studied) during the period in which the 

did not change. Initial rates in terms of (dA/dt),. were evaluated by 
plane mirror method while the pseudo first order rate constant k,, was 
calculated by using Guggenheim's method. The second order rate constants 
kj were obtained by dividing the k, by [S] wehre S is the reactant taken in 
excess. The pH of solutions was measured by using Systronics digital pH 
meter model-335. A high precision thermostatically controlled water bath 
with accuracy ± 0.1°C was employed for maintaining the constant 
temperature during the kinetic runs. The stoichiometry of the reaction was 
determined iodometrically. 

RESULTS AND DISCUSSION 

Reaction mixture containing 0.001 M [OCA] + 0.02 M [NalOJ + 2.5 
% (v/v) acetone, was prepared and kept for 24 hours. It developed violet 
colour changing into orange followed by precipitation. It was filtered after 
24 hours, filtrate extracted with petroleum ether and evaporated at room 
temperature. On separation of a droplet (due to water), a solid residue was 
left on the petri dish. This solid residue was dissolved in petroleum ether 
(40-60'’C) and a little amount of silica gel was added to it and mixed with a 
clean spatula so as to absorb the complete residue on silica gel. This silica 
gel was then applied on the top of a silica gel-H packed, 31 cm long column 
with 0.6 cm diameter. Gradient elution was employed using petroleum ether 
and chloroform in decreasing ratio. Four components were separated out of 
which a light orange colored major component was collected, checked for 
purity by using TLC and crystallized in ethyl alcohol. Its melting point was 
found to be 59® C. This compound was identified as a quinone [5]. The UV- 
VIS spectrum of this compound in CHClj, showed for this compound at 
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270 nm, 370 nm and 520 nm, which suggested the presence of quinonoid 
structure. 

Further, I.R. spectrum (in KBr) showed bands at 1632 cm '(s) (due to 
1,4-benzoquinone pattern); 3255 cm '(s) (due to overtone of C=0 strech); 
2713 cm '(w) (due to isolated C-H streching); 1514 cm '(s), 1469 cm '(s) 
and 1354 cm '(s) (due to C=C ring strech); 1200 cm '(s) and 1098 cm ' (due 
to in plane C-H bending); 825 cm‘'(m) (due to out to plane C-H bending in 
case of two adjacent H-atoms); 761 cm '(s) (a characteristic strong band 
due to C-Cl streching); and 612 cm '(m) and 518 cm '(m) both due to C-H 
bending mode). This product was finally charactrized as chloro-1,4- 
benzoquinone on the basis of its N.M.R. spectrum in CDClj which showed 
the peaks at 8 = 7.260 (s, IH) and at 8 = 6.991 (d, 2H) which are due to 
three protons of the ring [3,10,13]. 

Stoichiometry of the reaction was found to be I mol OCA : 2 mole 
periodate. The data presented in the table-1 were subjected to calculations 
by using the Vant hoff s differential equation and it was established that the 
reaction is second order,being first order in each reactant. The second order 
kinetics was also proved by the fact that the rate was linearly related to the 
concentration of the reactant varied in each case. Under pseudo first order 
conditions, the values of (dA/dt) ' vs [S] ' (Table-2) were also used to get a 
linear plot with almost negligible intercept, indicating that the intermediate 
formed in the slow step, should be unstable and it should have got consumed 
in a subsequent fast step. 

Table-1 _ 

^nuix = 460nm;pH = 5.0;*Acetone = 5.0%(v / v);?^ Acetone = 20%(v / v); 

Temp. = 30 ± 0.1" C 

[OCA] X 10’ M 50.050.0 50.0 50 . 0 ^ 50.050.0 ?t 3.0* 4.0* 5.0* 6.0* 7.0* 8.0* 

[NalOJx lO’M 5.0 6.0 7.0 8.0 9.0 10.0 30.0 30.0 30.0 30.0 30.0 30.0 

(dA/dtXxlO’ 0.70 0.85 0.95 1.00 1.25 1.375 3.3 3.8 4.8 6.0 7.0 8.0 

i 

(min') 
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Table-2 


= 460nni;pH = 

5.0;*Acetone = 

25.0%(v / v);^ Acetone 

! = 50%(v/v); 

Temp.= 

30±0.1°C 










OCAxKFM 

4.5* 5.2* 

5.4* 

5.6* 

5.8* 

6.0* 

0.3# 

0.3# 

0.3# 

0.3# 

0.3# 0.3# 

[NalOJx ItfM 

0.5 0.5 

0.5 

0.5 

0.5 

0.5 

3.0 

3.2 

3.4 

3.6 

3.8 4.0. 

(dA/dt).xl04 

3.4 3.5 

3.6 

3.8 

3.9 

4.0 

3.3 

3.5 

3.8 

4.0 

4.25 4.5 

(min-1) 












Kinetic studies were also made in the pH range 3.5 range 3.5 to 8.0 (Table- 
3) indicated that the rate was maximum at pH 5.0 The first part of the profile 
i.e. the increase in the rate from pH 3.5 to 5.0 may be to the decrease in the 
protonation of OCA from pH 3.5 to 5.0, which makes greater concentrations 
of OCA available for the reaction. This assumption is in lin with the fact 
that unprotonated OCA is the reactive species as shown by us in the 
mechanism proposed. Further, the concentration of the periodate monoanion 
is maximum around pH 5.0 and decreases beyond this pH value as worked 
out by earlier workers [4,6,11,15,16]. This may be the reason for the decrease 
in the rate of reaction beyond pH 5.0 This behavior also supports our 
assumption in the propsed mechanism (CHART) that out of the various 
species of the periodate, the species taking part in the reaction in present 
case is the periodate monoanion i.e. [lOJ. 


Table-3. Effect of pH on the reaction rate 


[OCA] = 3.0 X 10' M, [NalOJ = 3.0x10-^ M, = 460 nm Acetone = 5.0 % (v/v), 

Temp. 35.0 ±0.1 PC. 




pH 3.5 

4.0 4.5 5.0 5.5 

6.0 6.5 

7.0 7.5 8.0 

(dA/dt). X 104 20.0 

(min') 

25.0 30.0 40.0 30.0 

25.0 20.0 

15.0 10.0 5.0 
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The data in table - 4, suggests a linear relation between log (dA/dt)j and 1/ 
D (where D is the dielectric constant of the medium) with negative slope 
indicating an ion-dipole interaction in this reaction [1]. This point is well 
supported by the high negative value of the entropy of activation and our 
observation that the free radical scavengers like acry-amide and allyl alcohol 
exerted no effect on the reaction rate. A primary linear type plot between log 
(dA/dt)j vs ionic strength (p) was obtained by us by using the data given in the 
table-5, which also supports the ion-dipole interaction assumed by us in the 
mechanism proposed. 

Table - 4. Effect of dielectric constant of the medium on the reaction rate 

[OCA]=3.0x 10-'M,[NaIOJ=3.0x lO-'M, =460nm;pH=5.0,Temp.=30 ±0.1*^ 

Dielectric constant 72.4 71.0 69.7 68.4 

(dA/dt)i X103 (min ') 3.75 2.4 1.6 1.0 


Table - 5. Effect of dielectric constant of the medium on the reaction rate 

[OCA] = 3.0xia'M,[NaIO4] = 3.0x la^M, =460nm;pH = 5.0; 

Acetone = 5.0% (v/v); Temp. = 35.0 +_0.1°C. 

pxlO^ 4.0 5.0 6.0 7.0 

(dA/dt)i X loomin ') 3.5 4.0 4.5 5.25 


The values of different thermodynamic parameters evaluated from the linear 
Arrhenius plot obtained by using the results of the kinetic studies made at four 
temperatures ranging from 30 ± 0.1° C to 45 + 0.1° C and taking [OCA] = - 
0.003 M, [Nal OJ = 0.03 M and acetone = 10.0% (v/v) are, Ea = 10.94 k cal. 
mol '; A = 1.692 x 10° lit. mol ', sec-1 ; =- 32.11 e.u. ; ^*= 20.295 

k.cal.mol' and = 10.325 k.cal. mol'. Where, Ea is the energy of activation, 
A is the Arrhenius frequency factor, is the entropy of activation, is 
the free energy of activation and is the enthalpy of activation. 

A large negative value of entropy of activation is mainly observed in polar 
solvents and also suggests the formation of a charged and rigid transition state 
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which is expected to be strongly solvated in the polar solvent employed and is 
in accordance with the first and third steps in the mechanism proposed. The 
decrease in rate with the decrease in dielectric constant (Table-4) also supports 
this assumption. A low value of energy of activation is characteristic of a 
bimolecular reaction in the solution. 

On the bais of kinetic studies and the main reaction product characterized, 
the mechanism presented in the chart can be proposed which fits the rate law 
given below. 

(dA/dt) = kj [OCA][IOJ 



0 

Qilero-1,4 • b«aiMiiaDoaB 


mm-i\ [mi 
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The formation of a charged intermediate (I) which reacts with another 
periodate molecule to form quinoneimine (II) and subsequent fast hydrolysis 
of II to give chIoro-l,4-benzoquinone are the main features of the proposed 
mechanism is in good agreement with the reaction product characterized 
by us and the different kinetic results reported and discussed above viz. the 
order of reaction in each reactant, stoichiometry observed and the studies 
related to the effect of solvent, ionic, strength, pH, and free radical 
scavengers on the reaction rate as well as the thermodynamic parameters 
evaluated. 
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QUASI-CONTRACTIONS 
AND APPROXIMATE FIXED POINTS 
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(Received 20.12.2002) 

ABSTRACT 

In dtis note, we discuss approximate fixed points of quasi-contractions on metric spaces. 
Mathematics Subject aassifications (2000): 47H10, S4E40, S4H2S. 

Key wordsandphrases:Fixedpoint,iq)proxiniatefixed point,contiaction,quasi-contnictioa 


INTRODUCTION 

Let{X, r/) be a metric sp^ and 7? X-^X.Let ebeapositive number. Thenaz € 
A'is an 8 -fixed point of if d(Tz, z) < e • The mtp Tis said to have approximate fixed point 
property if, for each s >0> the map possesses at least one 8-fixed point (see, for 
instance, [10]). Thetheoiy of approximate fixed points finds applications in mathonatical 
economics, noncooperative game theory, dynamic programming, nonlinear analysis, 
variational calculus, theory of integro-differential equations and several other areas of 
j^plicable analysis (see, for instance, [ 1 ], [3], [4], [8], [10] and sev^ references thereof). 

The map TisaBanach contraction if 

d(Tx,Ty) < qd(x.y) for some 0 < ^ < 1 and allx,y inX (1) 

The map Tis called a quasi-contraction (cf Cfiric' [2]) if 

d{Tx, Ty) ^ q. max {dix, y), d{x. Tx), d(y. Ty), d(x, 2», d(y, Tx)) (2) 

for some 0^^ < 1 and all x, y vaX. 

We remark that T satisfying 

d{Tx, TyHq.m2X {d{x. y), d(x. Tx\ d(y. Ty), [d(x, Ty) + d(y. Tx]/2} (3) 

forsome0:g^<l andallx.yinA'neednotsatisfy(2).Noticethat(l) => (3) =>(2),and 
the reverse implication is not true. Indeed, among several generalizations of the Banach 
contraction, the quasi-contraction is considered the most general. For an excellent 
comparison of such maps, refer to Rhoades [9]. C'iric' [op. cit.] has shown that a 
quasi-contraction on a T-orbitally complete metric space has a unique fixed point. A 
metric space is called T-orbitally complete if every Cauchy sequence consisting of 
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elements of the Picard orbit 0(;Co,'»): ,..}of T for any is 

convergent to a point in X. It is well known that a quasi-contraction need not be 
continuous, while a Banach contraction is unifcmnly ccmtinuous. 

In all that follows q satisfying (1) or (2) will be called Banach contraction or quasi- 
contraction constant as the case may be. Recently Tijs et al. [10] have obtained the 
following af^roximate fixed point theorem for the Banach contraction. 


Theorem 1. Let (X, d) be a metric space and T a contraction on X with the Banach 
contraction constant q. Then, for each e > 0, FIX^ {T):{x € JT: d(7>f, jf) ^ g } # 0 and 

the diameter of FIX^ (7) is not larger than 2 e H\~q). 

The purpose of this note is to extend the above result to quasi-contractions and 
discuss its implications. 

APPROXIMATE nXED POINT THEOREMS 


The following is ourmain result. 

Theorem 2. Let (X, d)be a metric space and T a quasi-contraction on X with 
quasi-contraction constant q. Then for each e > 0, 

(I) FIXj {T)=={x 6 X:d(Tx, and (II) the diameter of FIX^ (f) is not larger 


than 


'ill 


\ 

e, 

y 


Proof. Pick Xg eX. Construct a sequence in X such that Tx^, n * 0,1,2,.... 

Then following C'iric'[2, p. 271], we have 


dirxg.rxg)< 




With m=/j +1 and n large, this implies that 


d{x„x„^^)< 




X,) <8. 


Therefore, for large n, x^ is an 8 -fixed point of T. This proves (I). 

Let w and z be any two 8 - fixed points of J, then 
diw,z)^ diw, Tw) + rf(7W. 7z) + diz. Tz) ^ 2 b+ d{Tw, Tz) 

< 28 + ^.max {d{yv, z), d{w, 7W), d{z. Tz), d(w, Tz), d(z,T\v)) 
^ 2z +q{d{w. z)+ 8}. 
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This yields (H). 

Corollary 3. Let the quasi-contraction Jin Theorem 2 possess a fixed point z. Then: 

(I) z is the unique fixed point of T; and 

(II) For each sequence {z^}withthepix>pertythat for eachne ^(naturals) the pointz^ is 
an 1 /n-fixed point, we have j[™ z^=z. 

Proof. (I) is immediate from (2). Since z is a fixed point, it must belong to fix (7) for 
each E > 0. Then, by Theorem 2 (II), ^ 

d(z^z) < diam (FIX„„(7))<(2+q)/{«(l.q)}. 

So, d{z^ z) = 0. This completes the proof. 

We remark that neith^ the complet^ess nor the orbital completeness of the space 
is required in all the above results. The following example illustrates our result. 

Example 4. Let X= {x: 0 < x < 1 andx is a rational numbo*} and Tx= 1/2, if x= 1 andO, 
otherwise. 

Then it is easily verified that T is a quasi-contraction with the fixed point 0, and T ^ (l)or 
(3) (see also Rhoades [9, Theorem 1 (xxvi)]). Further, for any 0 < s < 1/2, FIX (T) 
consists ofall points ofA"<^ [0,1/2). ® 
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